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ABSTRACT 

Context. The accretion history of low-mass young stars is not smooth but shows spikes of accretion that can last from months and 
years to decades and centuries. 

Aims. Observations of young stars in outbursts can help us understand the temporal evolution of accreting stars and the interplay 
between the accretion disk and the stellar magnetosphere. 

Methods. The young late-type star VI 118 Orionis was in outburst from 2005 to 2006. We followed the outburst with optical and 
near-infrared photometry; the X-ray emission was further probed with observations taken with XMM-Newton and Chandra during 
and after the outburst. In addition, we obtained mid-infrared photometry and spectroscopy with Spitzer at the peak of the outburst and 
in the post-outburst phase. 

Results. The spectral energy distribution of VI 118 Ori varied significantly over the course of the outburst. The optical flux showed 
the largest variations, most likely due to enhanced emission by a hot spot. The latter dominated the optical and near-infrared emission 
at the peak of the outburst, while the disk emission dominated in the mid-infrared. The emission silicate feature in VI 118 Ori is 
flat and does not vary in shape, but was slightly brighter at the peak of the outburst compared to the post-outburst spectrum. The 
X-ray flux correlated with the optical and infrared fluxes, indicating that accretion affected the magnetically active corona and the 
stellar magnetosphere. The thermal structure of the corona was variable with some indication of a cooling of the coronal temperature 
in the early phase of the outburst with a gradual return to normal values. Color-color diagrams in the optical and infrared showed 
variations during the outburst, with no obvious signature of reddening due to circumstellar matter Using Monte-Carlo realizations of 
star-l-disk-l-hotspot models to fit the spectral energy distributions in "quiescence" and at the peak of the outburst, we determined that 
the mass accretion rate varied from about 2.5 x 10~^ Mq yr~^ to 1.0 x 10~'' Mq yr~^; in addition the fractional area of the hotspot 
increased significantly as well. 

Conclusions. The multi- wavelength study of the VI 1 18 Ori outburst helped us to understand the variations in spectral energy distri- 
butions and demonstrated the interplay between the disk and the stellar magnetosphere in a young, strongly accreting star 

Key words, accretion, accretion disks - Infrared: stars - stars: circumstellar matter - stars: coronae - stars: pre-main-sequence - 
stars: individual (VI 1 18 Ori) - X-rays: stars 

1 . Introduction therefore, show a low recurrence rate. EXors (named after the 

prototype EX Lup), in contrast, show somewhat smaller out- 

The star formation process involves strong accretion of circum- bursts (A^ = 2-3 mag) on a much shorter timescale, from 

stellar matter onto the protostar. The time evolution of the mass ^ few mo nths to a few year s, and may occur repeatedly (see 

accretion rate is of deep interest to understand the timescale review by iHartmann & Ken von 1996; also Herbig 2008). Such 

of stellar growth and lifetime of proto-planetary disks. While outbursts are believed to originate during a rapid increase of the 

the mass accretion rate in young stars overall decreases with disk accretion rate over a short period of time, from values of 

increasing stellar age (Hartmann et al.||1998|), it can also show io-7 Mq yr^^ to IQ-^ Mq yr-\ although the underlying cause 

significant and rapid changes over time (e.g., [Hartmann et alj of such increase in mass accretion rate is unclear (thermal disk 

ll993|). A handful of accreting young stars display powerful erup- instabilities, Lin & Papaloizou 1985; Bell & Lin 1994; close 

tive events with flux increases in the optical regime of a few companions, Bonnell & Bastienlll992i ; cluster-in duced encoun- 

magnitudes. Two classes have emerged: FUors, which display ters.Pfalzner et al...2008l: gian t planets in the disk. lClarke &"S^ 



outbursts of 4 magnitudes and more, last several decades and, |1996l:lLodato & Clarkell2004l: a combination of gravitational in 



stability and the triggering of the magnetorotational instability. 

Send offprint requests to: e-mail: Marc . AudardSunige . ch Armitage et al. 2001; accretion of clumps in a gravitationafly 

* Tables lAjandO and Figure|lIT]are only available in electronic unstable disk. Vorobvov & Basull2005l 120061) . The limited num- 

format at http ://w w w. aanda. org , ber of eruptive young stars and the long recurrence time (espe- 



2 



M. Audard et al. : The 2005 Outburst of V 1 1 1 8 Ori 



cially for FUor-type objects) make it difficult to test models. It 
is, therefore, important to study in as much detail as possible the 
evolution of outbursts and to understand the place of outburst- 
ing young stars in the evolutionary scheme from highly accret- 
ing, embedded protostars to "normal" accreting classical T Tauri 
stars (CTTS). 

The optically revealed CTTS can make excellent compar- 
ison stars to outbursting young stars: they have been studied 
extensively in the past: optical and infrared observations have 
given clues about the presence of accretion disks, on the mass 
accretion rates, and disk winds, while millimeter observations 
provided constraints on the dust disk mass and on CO outflows. 
In the X-rays, the origin of emission in young, accr eting stars has 
gained significant attention in the last few years. iKastner et al.l 
(1200 2') obtained the first X-ray grating spectrum of TW Hya, a 
moderately accreting CTTS, and observed unexpected new spec- 
tral features: the O VII and Ne IX He-like triplets showed very 
low forbidden-to-intercombination line ratios and indicated high 
electron densities of the order of 10^^ — lO^'^ cm^''. The effect 
of UV photoexcitation on the observed line ratio was deemed 
negligible. In addition, the X-ray spectrum was consistent 
with a quasi-isothermal plasma of about 3 MK. iKastner et alj 
(I2002h concluded that the X-ray emission in TW Hya was due 
to acc retio n, a result furth e r su pported by IStelzer & SchmittI 
(12004') and ' Ness & SchmittI (12005). Additional high-resolution 
X-ray spectra of CTTS also indicated that accretion may play a 
significant role for X - rays in accreting star s (Schmitt et al. 20051 
Gtother e t all 120061; 'Robr ade & SchmittI 12006; .Argiroffi et all 
20071 iGunther et al. 20071: ISacco et al.1 l2008l) . Soft X-ray 



emiss i on from shocks in jets may a lso be detected (Giidel et al.' 
20051 120081; IKastne r et alj l2005t [R obrade & Schmitt 2007; 



Schneider & Schmitt 12008 : Gunther & Schmitt 2009) 



Telleschi et al., (.2007,) and iGiidel & TeUeschi. (,2007.) showed 



that CTTS display a soft X-ray excess from plasma at low 
temperature. Such a plasma (mostly detected in the O VII 
lines) is difficult to detect with CCD spectroscopy but can 
easily be revealed with high-resolution X-ray spectra when the 
absorbing column density is not too high. The soft X-ray excess 
can co-exist with the hot plasma observed in the vast majority 
of youn g stars, which is due to scaled-up solar-like magnetic 
activit y Ismifli et alJl2005yAudard et alj|2005al; iPreibisch et"!!] 
I200I iTeUeschi et alJ l2007r The origin of such an excess is 
unclear, but it could be due in part to accretion onto the stellar 
photosphere. It is crucial to understand how the X-ray emission 
in young stars, in particular those accreting matter actively, can 
be influenced by the accretion process. FUor and EXor stars are, 
therefore, the ideal cases to understand the physical mechanism 
common to both CTTS and outbursting stars. 

In the recent years, several young low-mass stars were 
observed in outburst, in the optical and infra red, but also 
in flie X-rays: VI 6 4 7 Ori i n 2003-2005 ( e.g.. IKastner et al J 
2004 lGrossoetal.1 I2005L IKastner et al] |2006at see also 



Aspin et al. 2008 and references therein) and recently in 2008- 
2009 (Itagaki 2008]; lAspiiil 120081; Fv^kat & Anandi^aol 120081: 
Aspin etal.ri2009l) . VI 1 18 Ori in 2005-2006 (lAudardetal 



2005bVLorenzetti et alj20"06ll2007l;lHerbig|2008l) . lKastner et al 



(2006b) also observed with Chandra a young stellar object in 
eruption in LDN 1415 but failed to detect i t. EX Lup also showed 
an extreme outburst in early 2008 (lJonesll2008a.b; Kospal et alj 
I2OO8I) . I Abraham et al.1 (l2009l) detected crystaUine feafiires in the 
silicate feature of EX Lup in outburst which were not present in 
the pre-outburst spectrum, suggesting crystallization by thermal 
annealing in the surface layer of the inner disk. X-ray obser- 
vations were obtained with Chandra (PI: Weintraub) and Swift 



(PI: Stringfellow) but are not yet published. We note also that the 
FUor objects FU Ori, V 1735 Cyg, and Z CMa were also detected 
in X-rays (Ski nner et a l. 2006, 2009a; Stelzeret al. 2009). 

The VI 647 Ori campaign showed an increase of the X-ray 
flux up to a factor of 200 from its pre-outburst flux, in line with 
the flux increase in the infrared. The X-ray flux then followed 
the optical outburst flux and returned to its pre-outburst level 
after the outburst ended (Kastneret al. 2006a). The initial ob- 
servations of V1118 Ori in X -rays showed a differ ent behavior, 
with little flux enhancement dAudard et al.ll2005bl) . but also in- 
dicated that the rapid increase of accretion rate in outbursts can 
impact the X-ray emission of young accreting stars. The present 
paper aims to present the remainder of the X-ray data of V1118 
Ori taken during and after the 2005-2006 outburst, together with 
contemporaneous optical and infrared data. 

2. The young erupting star V1 1 1 8 Ori 

The outburst of V1118 Ori, a young low-mass M le star in the 
Orion Nebula {d = 400 pc; iMuench et allllOOSl for a discus- 
sion), was reported i n earl y Jan uary 2005 by Williams et al.l 
( 120051) . iHiUenbrandl d 19971) and IStassun et all (ll999|) provide 
details about its physical properties: M^, = 0.41Mq; Ri, = 
1.29i?Q; Prot = 2.23 ± 0.04 d; Lboi > 0.25 Lp,;\osT^s\K ] = 
3.562; log i[?/r] = 6.28. Recentlv. iReipurth et al.r(l2007l) re- 
solved V1118 Ori into a close binary separated by O'.'IS with 
a position angle of 329° and a magnitude difference of Am — 
0.4 mag in the Ha-band image. The observation was obtained 
on 2004 Jan 29, thus before the 2005-2006 outburst (see Fig.|6ll. 
The binarity of V1118 Ori leaves unclear which star actually 
erupted. However, we note that the small magnitude difference 
indicates that the components are similar, suggesting similar ef- 
fective temperatures. We may also assume that they show similar 
disk evolution. While the binarity of VI 118 Ori complicates the 
interpretation of the combined photometry and spectroscopy in 
quiescence, this should have little impact on the interpretation 
of the outburst data, as only one star-ndisk component dominates 
the emission. In this paper, we have used the above stellar prop- 
erties that assumed a single star 

VI 118 Ori has shown freq uent outbursts (e.g., 198 3-84, 
1988-90, 199 2-94, 1997-98; see iGarcia & Parsamianl |2000 and 
lHerbigll2008l for details). In fact, after returning in quiescence 
in mid-2006, VI 118 Ori had another outburst in late 2007 
(iGarcia & Parsamianll2008l) . We focus here on our monitoring 
campaign to study the 2005-2006 outburst of VI 1 18 Ori in the 
X-rays, optical and infrared. Additional properties in the opti- 
cal and infrare d were obtained independen tly d uring and afte r 
the outburst bv lLorenze tti et al. (2006, 2007) and Herbig (2008!). 
Using wind models. lUorenzetti et al. (2006) derived a mass loss 
rate of 4 x 10~* Mq yr~^ from the H I recombination line, and 
(3 - 8) X 10"'^ Mq yr~i from the CO emission at 2.3 /xm for the 
neutral molecular gas; they also found n o evidence of infrared 
cooling from a collimated jet or outflow. ILorenzetti et al.l (|2007|) 
also found evidence of intrinsic polariza tion in the / ban d, and ir- 
regular fluctuations during the outburst. iHerbid (l2008l) obtained 
Keck/HIRES spectra of VI 118 Ori in the decaying phase of 
the outburst and in the post-outburst phase. He noted the de- 
tection of Li I A6707 in emission during the outburst, in contrast 
with the absorbed feature in CTTS spectra. The feature, how- 
ever, returned in absorption after the outburst. A similar behav- 
ior occurred in the K I lines at AA 7664 and 7698 A. A P Cyg 
profile was also found in Ha during th e outburst, which disap- 
peared thereafter In the near-infrared, ILorenzetti et all (120071) 
also found that the emission lines detected during the outburst 
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Table 1. X-ray observation log. 



Parameter Sep 2002 Jan 2005 Feb 2005 Mar 2005 Sep 2005 



Satellite Chandra Chandra XMM-Newton XMM-Newton XMM-Newton 

ObsID 2548 6204 0212480301 0212480401 0212481101 

Duration (ks) 48 5 20 20 20 

Observation date 2002 Sep 6-7 2005 Jan 26 2005 Feb 18-19 2005 Mar 21 2005 Sep 8 

UT 12:57-02:54 03:07-05:04 22:36-04:50 16:20-22:08 12:43-18:51 

Average JD - 2,450,000 2,524.3 3,396.7 3,420.6 3,451.3 3,622.2 



Jan 2006 Feb 2006 Mar 2006 Apr 2006 Jul 2006 Dec 2007 

Satellite Chandra Chandra XMM-Newton Chandra Chandra Chandra 

ObsID 6416 6417 0403200101 6418 6419 8936 

Duration (ks) 30 30 100 30 30 37.5 

Observation date 2006 Jan 4-5 2006 Feb 23 2006 Mar 2-3 2006 Apr 23 2006 Jul 24 2007 Dec 14-15 

UT 18:33-03:09 09:41-17:47 19:02-21:06 14:42-23:31 07:50-16:49 14:49-01:13 

Average JD - 2,450,000 3,740.5 3,790.2 3,797.8 3,849.3 3,941.0 4,449.3 



(Hi, Hex, CO 2.3 /im band, a few neutral metals) disappeared 
about a year later. 

In the X-ray regime, the i nitial Jan-Mar 2005 data were pub- 
lished bv lAudard et"ai] (l2005bi) . In brief, the X -ray data of early 
2005 indicated that the X-ray flux and luminosity stayed simi- 
lar within a factor of two during the outburst, and at the same 
level as in a pre-outburst observation in 2002. The fluxes in 
the optical and near-infrared varied more significantly, within 
factors of 2 — 10. The hydrogen column density showed no 
evidence for variation from its modest pre-outburst value of 
N}i « 3 X 10^^ cm^^. However, there was evidence of a spec- 
tral change from a dominant hot plasma (« 25 MK) in 2002 
and in January 2005 to a cooler plasma (« 8 MK) in February 
2005 and probably in March 2005. We hypothesized that the hot 
magnetic loops high in the corona were disrupted by the clos- 
ing in of the accretion disk due to the increased accretion rate 
during the outburst, whereas the lower cooler loops were proba- 
bly less affected and became the dominant coronal component 
dAudard et al.l l2005bl) . We argued that the cool component in 
VI 1 18 Ori could not originate from shocks because free-fall ve- 
loci ties of matter falling from the truncation radius are too low 
(see lAudard et al.ir2005b. for f urther discussion). In a subsequent 
paper, iLorenzetti et al. 1 1200^ ) independently analyzed our pub- 
lic data sets of early 2005, and also included the September 2005 
XMM-Newton observation, together with their near-infrared data 
sets. The September 2005 observation showed a decrease in X- 
ray flux at the start of t he decay phas e. 

The initial paper bv lAudard et al.l (l2005bl) presented and ana- 
lyzed the outburst data through March 2005. We present here the 
XMM-Newton and Chandra data obtained from September 2005 
on, which are analyzed in the context of our multiwavelength 
analysis. We also include a post-outburst data set obtained by 
Chandra in December 2007. The latt er was taken during a mi - 
nor outburst detected in the optical (Garcia & Parsamianll2008l) . 
Extensive optical and near-infrared photometry are presented 
and analyzed, obtained from our team's imaging data, along with 
mid-infrared photometry and spectroscopy with Spitzer. 



3. Observations and data reduction 

3.1. X-ray 

Table [Uprovides the observation log of our 2005-2006 monitor- 
ing campaign of th e outburst of VI 1 18 Ori with XMM-Ne wton 
dJansen et al.ll2001l) and Chandra jWeisskopfet al .1119961) . We 



also provide the information about the 2002 serendipitous obser- 
vation of VI 1 18 Ori (an obser vation i n 2001 with XMM-Newton 
was reported in I Audard et aLll2"005bl but is not mentioned here 
since the star was not detected). A deep XMM-Newton obser- 
vation was obtained in early 2006 that complements the short 
monitoring observations. Finally, we include the post-outburst 
Chandra observation in December 2007 as well. Note that the 
angular resolutions of both X-ray satellites were not high enough 
to separate the V 1 1 1 8 Ori binary. 

The XMM-Newton data were processed with SAS 7.0. 
Standard procedures were applied. We used an extraction cir- 
cle of radius 20" for the source and a nearby background cir- 
cular region of 60" radius (40" for Sep 2005 and 35" for Mar 
2006). Event patterns lower than 4 and 12 were used only for 
the European Photon Imaging Came ras (EPIC) pn and MOS 
(IStruder et al.ll200ll: iTumer et al.l l2001 ), respectively. The back- 
ground flux levels were high during all XMM-Newton observa- 
tions, and in particular i n Marc h and September 2005. As de- 
scribed in lAudard et"aLl (l2005bl) . we used only the MOS 1 and 
MOS2 data for the March 2005 observation. The September 
2005 observation was so affected by the background that we lost 
about half the exposure time in the EPIC pn, while there was 
no MOS data available. The deep March 2006 observation was 
affected by a system failure at the Mission Operations Centre 
and the EPIC pn experienced full scientific buffer in the last part 
of the observation, explaining the reduced exposure time in the 
EPIC pn (78ks) compared to the EPIC MOS (90 ks). No Optical 
Monitor data were taken with the XMM-Newton X-ray observa- 
tion due to the presence of the nearby bright Trapezium stars. 

The Chandra data were processed with CIAO 3.3 and 
CALDB 3.2.3B The task psextract was used to extract the spec- 
tra for VI 1 18 Ori and the nearby background. We used a circle 
of 2" radius (« 4 pixels) for the star and an annulus centered 
at the position of the star but with radii of 10 and 60 pixels for 
the background; our background area was, therefore, 212 times 
larger than the source extraction area. Our extraction radius for 
the star includes 95% of the encircled energy at 1.5 keV and 
90% at 4.5 keV. Note that for the January 2005 observation in 
sub-array mode, we used an outer radius of 40 pixels for the 
background annulus. For the 2002 observation, we used two cir- 
cles of radii of 1 5" and 45" for the so urce and the background, 
respectively (see lAudard et al.ll2005bl) . 



We used the pipeline data for the December 2007 observation, 
which was calibrated with CALDB 3.4.2. 
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Table 2. Spitzer IRAC and MIPS flux density measurements (units of mJy). 



Date 3.6 /im 4.5 /im 5.8 fim 8.0 fim 24 /im 



2004 Mar 09 28.6 ±0.1 36.3 ±0.1 39.2 ± 1.3 54.8 ± 0.5 

2004 Mar 20 ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ 74.2 ± 0.5 

2004 Oct 12 46.5 ±0.2 50.6 ± 0.2 53.3 ± 0.3 55.9 ±0.5 

2004 Oct 27 •■■ 43.4 ±0.1 •■■ 57.5 ± 0.3 

2005 Feb 20 143.3 ± 1.7 169.7 ± 2.0 148.2 ± 4.9 144.1 ± 1.9 

2005 Mar 28 173.7 ± 1.8 203.9 ± 2.0 194.1 ±5.0 167.6 ± 1.9 

2007 Oct 21 42.4 ± 1.7 48.7 ± 1.9 52.2 ± 4.9 53.1 ± 1.8 



3.2. Optical and near-infrared 

3.2.1. SMARTS 

We observed V1118 Ori with the ANDICAM dual-channel im- 
ager on the SMARTS/CTICEI 1 .3m telescope on 115 nights. Such 
dense coverage is possible by the service mode operations em- 
ployed by the SMARTS consortium. We have 5 pre-outburst ob- 
servations (2004 Feb 2 through 2004 Apr 8) that can be used 
to set the quiescent flux levels. The intensive monitoring began 
2005 Jan 10, shortly after the outburst was reported. There were 
no observations from 2005 Apr 4 through 2005 Jul 31 due to 
proximity to the Sun. We observed with cadences between 1 
per day and 1 per week, with observations generally every 3- 
4 days. Our intensive monitoring terminated on 2006 May 6. We 
also followed the decay of the second outburst from 2007 Dec 
27 through 2008 Mar 27. Details about the data reduction are 
given in AppendixlAl Table lA.2l (available online only) gives the 
nightly average photometry for SMARTS. 

3.2.2. Villanova 

Photometric coverage in the optical (standard Bessel VRI) was 
obtained with the Celestron 14" optical tube assembly with 
a Paramount ME German equatorial mount at the Villanova 
University Observatory located on the campus near Philadelphia, 
PA. Observations were carried out with a SBIG ST7-XME de- 
tector thermoelectrically cooled to -25 degrees Celsius. Dark and 
flat field frames were collected at the end of each nights obser- 
vations. A standard error of 0.05 mag was estimated from the 
signal-to-noise ratio and seeing conditions. Table |A3] (available 
online only) lists the magnitudes obtained at Villanova. 

3.2.3. Additional data 

We have used published optical and near-i nfrared photometric 
data from lLorenzetti et al I J20q7l) (J JgJiD. lGarcia et al.1 ( l2006l) 



and lGarcia & Parsamianl (l2008lui^V 



3.3. Spitzer 

In addition to our observations of V 1 1 1 8 Ori in outburst (pro- 
gram ID 3716, PI: G. Stringfellow) and in post-outburst (pro- 
gram ID 41019, PI: M. Audard), VI 1 18 Ori was ser endipitously 
obser ved with the InfraRed Array Came ra (IRAC: [ Fazio et al] 
12004 by the Spitzer Space Telescope ( Werner et all |2004|) in 
March 2004 and twice in October 2004 (program IDs 43 and 
50, PI: G. Fazio). We show the IRAC images taken before the 
outburst in March 2004 centered on VI 118 Ori (Fig. [TJ. We 



SMARTS, the Small and Medium Aperture Research Telescope 
Facility, is a consortium of universities and research institutions that 
operate the small telescopes at Cerro Tololo under contract with AURA. 



prov ide the IRAC fluxes f or all observations in Table |2] (see 
also iLorenzetti et al.1 l2007l for the IRAC data of October 27, 
2004). Details about the Spitzer IRAC data reduction are given 
in Appendix IB] 

VI 118 Ori was also observe d with the M ultiband Imaging 
Photometer for Spitzer (MIPS; iRieke et al.l 12004) at 24 ^m 
before the outburst, on March 20, 2004 (program ID 58, PI: 
G. Rieke). Again, Appendix IB] provides the details of the data 
reduction. No MIPS photometry is available for VI 1 18 Ori dur- 
ing or after the outburst. Our programs (3716 and 41019) in- 
cluded MIPS spectral energy distribution data {R « 20) in the 
70 /im band; however, the on-time exposure (180 s) did not al- 
low us to detect VI 1 18 Ori, even during the outburst. Indeed, the 
background level (due to diffuse emission in the Orion nebula) 
produced a much higher flux (of order 40 Jy at 70 /im) than the 
expected signal from VI 1 18 Ori (of the order 0.1 Jy if extrapo- 
lating the SED at 70 /tm, see Fig.[T2b. 

Spitzer also observed VI 118 Ori with the InfraRed 
Spectrograph (IRS; Houck et alj|2004l) twice during the outburst 
(PID 3716) with the Short-Low (SL: 5.2-14.7 /tm, R ^ A/AA « 
64), Short-High (SH: 9.9-19.6 R = A/AA « 600), and 
Long-High (LH: 18.7-37.2 /im, R = A/AA « 600) modules 



IRAC1 



IRAC2 



IRAC3 




Fig. 1. Spitzer IRAC images centered on VI 1 18 Ori, taken pre- 
outburst in March 2004. The images are shown with a linear 
scale from to 200 MJy sr"^ The Herbig Ae star V372 Ori is 
seen near VI 1 18 Ori. 
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Fig. 2. Spitzer IRS background-subtracted spectra obtained near the peak of the outburst (panels a and b) and after the outburst 
(panels c and d). Panel d shows the SL and LL module data. Panels a, b, and c show the SL spectra as thick black lines, whereas the 
SH and LH spectra are shown as red and blue hnes. Detected emission lines are labeled (while other lines well-subtracted by the 
background spectrum are shown in italics). 

Table 3. Tentative Spitzer IRS line fluxes derived from the SH and LH spectra. 



Line 



A(^m) Flux(10"^** Wm^ 

Feb 2005 Mar 2005 



Nov 2008 



Ha 


- 5'(2) 

'Pl/2-'P3/2 


12.2786 






7.8 ±0.8 


[Neil 


12.8136 


41.1 ±0.7 


39.7 ± 1.1 


42.7 ±0.9 


[Sin] 


3P2 - 'Pi ... . 


18.7130 


46.5 ± 1.4 


42.5 ± 1.5 


34.6 ± 1.2 


H2O 


- 5'(0) 


28.2188 


15.2 ± 2.2 


18.9 ±3.4 


9.5 ±1.6 


[Siii] 


'Pi-'Po.... 


33.4810 


463.3 ±9.6 


374.4 ± 7.9 


183.4 ±4.7 



(2005 February 18 and 2005 March 11). No background ob- 
servations were taken with the high-resolution modules. On the 
other hand, the post-outburst data (PID 41019; 2008 November 
14) were taken with background spectra for the high-resolution 
modules and also included the Long-Low (LL: 14.0-38.0 /im, 
R = A/ A A w 64) module. 

Figure |2] shows all IRS spectra after background subtrac- 
tion. We describe in Appendix |B] our methodology to derive 
the background-subtracted spectra of VI 118 Ori. We also dis- 
cuss in the Appendix the validity to use the post-outburst back- 
ground observation for our outburst IRS SH spectra. In brief, 
the post-outburst, background-subtracted high-resolution spec- 



trum is well-subtracted for the continuum, as the flux is con- 
sistent with the low-resolution spectrum. In the case of the out- 
burst SH spectra, the continuum is also accurate below 14 /im 
since they are consistent with the low-resolution SL spectra. 
However, we prefer to stay on the safe side and claim that the 
observed increase in continuum flux for A > 14 fim during 
the outburst is unreliable, since we have no MIPS photometry 
or low-resolution spectra to confirm the increase. We also urge 
caution with regards to the detection of lines in the background- 
subtracted SH and LH spectra. The strong background line emis- 
sion of the Orion nebula, and its inhomogeneity (see Fig. IB. 2b 
make the background subtraction difficult, although not impos- 
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Fig. 3. Zoomed in regions of the Spitzer IRS background-subtracted spectra. The post-outburst spectra are shown in the bottom 
panel, while the top two panels are for the outburst spectra. In the leftmost panels, we show both the high-resolution SH (thin, 
red) and the low-resolution SL (thick, black) data. Detected emission lines are labeled in the bottom panels (while other lines well- 
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sible (e.g., a very strong PAH emission at 11.3 yum is well- 
subtracted). The[S III] emission lines at 18.71 and 33.48 iim are 
likely robust detections as their detection is consistent with the 
low-resolution SL and LL data. The [Si ll] is also generally well 
subtracted. Interestingly, the rotational H2 S(0) line is detected 
during the outburst and in post-outburst. However the S(l) line is 
not detected (it is is even slightly oversubtracted in the outburst 
spectra), and the S(2) line is only detected in the post-outburst 
spectrum (a faint excess in the outburst spectra is present only in 
one bin, likely due to slight differences in the wavelength scales 
of the source and background spectra). It is likely that the H2 
emitting conditions vary spatially in the Orion nebula, and that 
the background spectrum did not accurately reflect the H2 back- 
ground line fluxes near VI 118 Ori, creating an excess in line 
emission observed in the post-outburst spectrum. Alternatively, 
if the H2 emission lines are real, at least in the post-outburst 
spectra, the physical conditions near V 111 8 Ori are such that 
the para (odd quantum number J) lines dominate the emission. 
A similar case of difficult background line subtraction arises in 
the [Ne ll] line at 12.81 fim: its detection in the SH spectrum is 
also subject to caution, as it is not detected in the low-resolution 
spectra (although it should have been, see the discussion in the 
Appendix). In any case, we provide in Table[3]the tentative, mea- 
sured line fluxes from the high-resolution module spectra, while 
Figure [3] shows zoomed in regions of the IRS spectra of the de- 
tected lines. 

4. X-ray analysis 

4.1. Light curves 

Figure |4]presents the X-ray light curves for the Chandra obser- 
vations. The lower panels show the binned light curves, whereas 



the upper panels show the photon CCD energy as a function 
of time. Except for Sep 2002, the background is not subtracted 
in the light curves and not shown in the upper panels, since it 
contributed to less than 1 count to the Chandra total (source + 
background) number of counts over the total observing duration. 
Since VI 1 18 Ori was observed serendipitously in the Sep 2002 
observation and was placed at the edge of the ACIS-I camera, the 
background contributed more significantly; therefore, the light 
curve shown in the bottom panel is background-subtracted, and 
the CCD energy vs time background events are shown together 
with the CCD energy vs. time total events (for clarity, every 
ninth background event was plotted, which corresponds approx- 
imately to the area-scaled background contribution). 

Small-scale variability is observed in the Chandra light 
curves; however, no strong flares were detected during the ob- 
servations. The January 2006 observation displayed nevertheless 
a significant increase in the count rate in the last 10 ks, probably 
due to a moderate flare. The average count rates are 3.14, 2.31, 
1.10, 0.39, 1.23, 1.23, and 1.09 ct ks^^ (0.2 - 4.0 keV, except 
for Sep 2002: 0.2 - 6.0 keV) for the Sep 2002, Jan 2005, Jan-Jul 
2006, and Dec 2007 observations, respectively. While the effec- 
tive areas in all pointed observations are similar, the Sep 2002 
effective area was about 40-50% lower (due to the off-axis posi- 
tion of VI 1 18 Ori). Therefore, to compare the above Sep 2002 
count rate with the later Chandra count rates, one needs to mul- 
tiply by a factor of about 2, i.e., 6.00 ct ks^^. Clearly, assuming 
only a change in emission measure, the X-ray flux dropped by 
a factor of 4-6, and even 15 in Feb 2006, during the 2005 out- 
burst compared to 2002. It is important to emphasize that the 
serendipitous observation did not show evidence of strong flar- 
ing, and VI 1 18 Ori was likely caught in "quiescence". 

Figure |5] presents the X-ray light curves for the XMM- 
Newton observations. In comparison with the Chandra data, the 
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Fig. 4. Light curves for the Chandra observations (0.4-6.0 keV 
range). The upper panels show the CCD energy of the events as 
a function of time after the start of the observations, whereas the 
lower panels show count rate light curves (fj,) with a bin size of 1 
ksec. The average count rates are also shown as horizontal dotted 
lines. The time span for each panel was kept similar and equal 
to the time span of the Sep 2002 observation 50 ksec). Since 
the other observations were shorter than 50 ksec, the rest of the 
time span is marked as hashed regions. See text for a detailed 
description. 



XMM-Newton data were heavily impacted by background radia- 
tion. The (scaled) background contributions generally were sim- 
ilar or even higher than the net source contribution. In particular, 
during the March 2005 observation, the background completely 
overwhelmed the EPIC pn data which could not be used. No 
obvious large flare was observed during the XMM-Newton ob- 
servations, like in the Chandra observations. However, in March 
2006, VI 118 Ori showed a significant increase in X-ray flux in 
the last 35^0 ks of the observation, probably due to flares or the 
onset of an active region. This increase is observed in all EPIC 
cameras, comforting us that this behavior is not due to an in- 
crease in background flux that was improperly subtracted in the 
EPIC pn data. 

The light curves in Figures |4] and |5] suggest that the spectral 
fits represent snapshots of the thermal emission measure distri- 
bution of V1118 Ori's corona observed serendipitously before 
the 2005 outburst, and over the course of the optical/infrared 
outburst. 
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Fig. 5. Background-subtracted count rate light curves (/i) for 
the XMM-Newton observations, with bin sizes of 500 s and in 
the 0.4-6.0 keV range. The EPIC pn data are shown, except in 
March 2005 for which we summed the data of both MOS detec- 
tors. The (scaled) background hght curves are shown as dashed 
curves. 



4.2. X-ray and optical/near-infrared flux correlations 

Figure |6] shows the long-term light curve of V1118 Ori in the 
X-rays, optical, and infrared, from its pre-outburst detection in 
September 2002, through its outburst in 2005 and 2006 to the 
post-outburst phase at the end of December 2007. We have 
looked into correlations between the X-ray and optical and in- 
frared photometry. Since we obtained a few sequential optical 
and infrared photometry points every night of our campaign, it 
was not possible to determine short-term variations to compare 
with the X-ray light curves. Nevertheless, to quantify further the 
correlation between the optical and near-infrared flux densities 
with the X-ray flux, we have calculated average optical and near- 
infrared flux densities within ±20 days of the 6 X-ray observa- 
tions between January 2005 and February 2006 (we did not in- 
clude the March 2006 observation, as this one might have been 
contaminated by a flare). The optical and near-infrared flux den- 
sities are well-correlated with the X-ray fluxes measured at Earth 
(Fig- 13 as indicated by Pearson's correlation coefficients, which 
are in the range of p = 0.8 — 0.95. The sparsity of X-ray data 
points does not allow to measure any clear delay between the X- 
ray and optical/near-infrared flux densities. We have also done 
cross-correlation studies with the optical and near-infrared data 
(including only the data points for which we had simultaneous 
measurements). The onset of the outburst is unclear in most light 
curves (except perhaps in the V band), but the end of the outburst 
occurs at a similar epoch for all bands (around MJD 53860). 
The shape of the light curves also differ: the fluxes vary most at 
short wavelengths (BVRI), while the near-infrared fluxes dis- 
play shallower variations and they reach their peak some time 
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later (e.g., K peaks around MJD 53470 while V peaks around 
MJD 53433). The overall shapes of the light curves suggest that 
the mechanism dominating in the optical and near-infrared bands 
may be the same (e.g., due to a hot spot), although the different 
peak times suggest contamination of another mechanism long- 
ward of 1 ^m, probably disk thermal emission (see below). 



the me dian quantile Qso, and the index y — 2> x Q25 / Q75 in the 
t;-axis. lHong et al.l (|2004|) define the a% quantile, Q^,, as 



4.3. Spectral fits 

lAudard et all (12005 ^) presented spectral fits of the V1118 Ori 
data taken through March 2005. Lorenzetti et al. (2006) repeated 
the analysis for the February 2005 data and added the September 
2005 data, noting a decrease in the X-ray flux of V 1 1 1 8 Ori at 
that period compared to earlier in the outburst. For this paper, 
we have reprocessed the 2005 data with the latest calibration 
(see Sect. O. 

We used XSPEC 11 (jAmaud 1996) to fit the Chandra 
ACIS and the XMM-Newton EPIC spectra for each epoch. 
For February 2005 and March 2006, we fitted the EPIC pn, 
MOSl, and MOS2 spectra simultaneously, while we used the 
EPIC MOS only for March 2005, and the EPIC pn only for 
September 2005 (no MOS data were available). In general, we 
used a l-T collisi onal ionization equilibrium (CIE) model (apec; 
ISmith et dll2001h with a photoelectric absorption model, except 
for the deep March 2006 observation for which the signal-to- 
noise was large enough to use a 2-T CIE model. We have fitted 
the plasma metallicity in the high signal-to-noise ratio spectra 
of February 2005 and used the best-fit value, Z = 0.17Zq, for 
the other epochs. The coronal a bundances are relative to th e so- 
lar photospheric standard set of Grevesse & Sauval (Il998i) . The 
best-fit metallicity is in line with the values meas ured for Fe in 
the co ronae of young stars in Orion (Table 3 in iMaggio et alj 
12007 1). Since the spectral fits of the February and March 2005 
spectra are similar to those reported in Audard et al. (2005b), 
we provide in this paper their best-fit values (but adapt the emis- 
sion measure and luminosities to the adopted distance of 400 pc). 
Note that the September 2005 data are heavily contaminated by 
the high back ground level during the observation. Contrary to 
iLorenzetti et a l. (2006), we preferred not to provide spectral fits 
for this observation; however, we provide estimates for the ob- 
served X-ray flux at Earth and absorption-corrected X-ray lu- 
minosity, based on the previous February 2005 and the posterior 
March 2006 XMM-Newton observations. The results of our spec- 
tral fits are given in Table |4] 

4.4. Quantiles 

In view of the low observed count rates with XMM-Newton 
and Chandra during the outburst, we have explored an alter- 
native method to determine spectral properties for low-count 
X- ray spectra and have used the quantile analysis presented 
by iHong et al] (|2004J^ . This type of analysi s was, in particu- 
lar, successfully used by iGrosso et~a!] (l2005h for their analysis 
of the spectral properties of VI 647 Ori (see also ^ Skinner et alj 
l2009bl) . In brief, this method makes direct use of the event en- 
ergy values, and determines spectral properties based on quan- 
tiles of the total number of counts (e.g., median 50% quantile, 
and quartiles, 25% and 75%), instead of spectral energy bins. 
Such quantiles are then used as indicators for the X-ray color of 
the source. In particular, we make use of a diagram with the x = 
\ogiQ{m/ (1 — m)) index in the x-axis, where m corresponds to 

The code is available at |http : //hea-www . harvard. edu/C h, 
We used version 1.7. 



where E^p and £"10 are the lower and upper boundaries of 
the used energy band, and Ea% is the energy below which the 
net counts is a% of the total number of counts. Note that this 
method is able to take the distribution of background event en- 
ergies into account. We refer the reader to lHong et al.l(l2004l) for 
more details. 

The total (source + background) and background extraction 
regions used in the quantile analysis were the same as those for 
spectral fits. The energy bands used to determine the extraction 
of events are reported in Table|4]together with the 25%, 50%, and 
75% quantiles, the EgQ% energy, and the x and y indices. We also 
provide estimated values for the hydrogen column density and 
plasma temperature derived from the x and y indices. Finally, 
Figure[8]shows the color-color diagrams for all observations. 

There exist s a certain degree of degeneracy in the x and 
y indices (e.g.. iGrosso et al.ll2005h : indeed, the same color in- 
dices can sometimes be described by different combinations 
of column density A'h and plasma temperature T. However, 
these degeneracies are limited to certain sets of T and A^h: for 
example, in the September 2002 observation, color indices of 
X ~ —0.8 and y = 1.8 can represent any plasma tempera- 
ture from about kT w 0.2 to about kT « 0.8 keV and A^h 
can range from to 0.75 x lO^'^ cm^^. At higher temperatures, 
the degeneracy breaks down, and the observed x/y values of 
VI 1 18 Ori can be better constrained to fcT « 1.2-2.5 keV and 
A^H ~ (1 — 8) X 10^^ cm~^. The degeneracies differ slightly 
depending on the response matrix generated in the observation. 
For example, with our Chandra monitoring data, the degenera- 
cies occur around {x,y) ~ (—0.7,1.85) and at (—0.4,1.65). 
The first set again indicate a degeneracy for kT < 0.8 keV and 
A^H < 1 X 10^^ cm^^. The second set shows that the plasma 
properties cannot be disentangled between 0.8 and 1.5 keV for a 
column density of (3 — 8) x 10^^ cm^^. 

In contrast to the temperatures obtained from spectral fits, 
the temperature derived from the color indices are more sta- 
ble and vary little from 10 — 20 MK, except in February 2006 
{T^y « 35 MK) and March 2006 {T^y w 90 MK). In the case 
of March 2006, the quantile temperature is in fact close to the 
second temperature component (which also has the higher emis- 
sion measure) obtained with the 2-T fit. In the case of February 
2006, the low number of counts did not allow a good constraint 
on the temperature. There also appears to have some kind of 
degeneracy: while the spectral fit gives a low plasma tempera- 
ture and "high" column density, the quantile analysis gives the 
inverse. It remains unclear which result is more robust, but the 
quantile diagram does not show evidence of a degeneracy. It is 
then likely that the spectral fit reached a local minimum at a low 
temperature and could not find a combination of high plasma 
temperature and low A'h. The February 2005 quantile analysis 
and spectral fit both confirm that the temperature of the coro- 
nal plasma was cool er than in Sep t ember 2002 and March 2005, 
as initially report in lAudard et al.l (l2005bl) . The error bars of the 
color indices in January 2005 do not allow us to constrain the 
temperature of the coronal plasma, while the spectral fit with 
iiCttfitS/(P(ifi8tes^io evidence of cool plasma at the start of the 
outburst. 
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4.5. X-ray parameter correlations 



We have attempted to search for correlations between X-ray 
properties (iVn, T, N^^^y, T^y, EM, Lx, £'5o%)' but found gen- 
erally little or no correlation based on the Pearson correlation 
coefficient, p. Indeed, its absolute value amounted to values typ- 
ically below 0.5. One exception was a mild correlation between 
the average median photon energy, E^q%, and the absorbing col- 
umn density derived from the quantile analysis, N-^^xy, which 
showed p — 0.57. This is expected, as larger column densities 
will absorb low-energy photons. However, we emphasize that we 
have only used the estimated Nn^xy based on the location of the 
xy values in the quantile color-color diagrams, and that we have 
not taken into account the (relatively) large range of column den- 
sities covered by the quantile errors (typically, 1 — 5x 10^^ cm^'^, 
similar to values derived from spectral fitting). Another excep- 



tion was a strong correlation between the emission measure and 
the X-ray luminosity, with p = 0.96. Together with the above 
lack of correlation between the plasma temperature and Lx, and 
A^H and Lx, this is a strong indication that the X-ray variabil- 
ity observed during the outburst of VI 1 18 Ori is only due to the 
amount of material in the corona. It is, thus, probably related to 
the amount of mass falling from the disk into the stellar magne- 
tosphere, i.e., to the mass accretion rate. 

5. Optical and infrared analysis 

5.1. Silicate feature 

Figure |9] shows the silicate feature obtained at the three different 
epochs and its normalized version (Si, = 1 + F^/Fc). We have 
use a linear fit to the SL spectra in the log A — log plane, using 
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only wavelength ranges from 6.0-8.0 /im and 13.0-13.5 fim to 
determine the underlying continuum flux level, Fc- The shapes 
of the silicate feature are compatible for all three observations, 
but the Mar 2005 silicate feature appears to be brighter than the 
Feb 2005 or post-outburst (Nov 2008) normalized fluxes, i.e., 
more flux coming from the less optically thick disk upper layers, 
likely due to hotter temperature of such layers. The heating of the 
upper layers is probably due to the strong irradiation, e.g., by a 
hot spot (see below). The ratio S'n.a/S'g.s is equal to 0.95 — 1.0, 
suggesting grain growth, and the ratio (S'n.s/S'g.s)/ (S'p"^^'") is 
0.74 — 0.79, evidence of a flat silicate feature. 

In the outburst spectrum of the EXor prototype EX Lup, 
[Abraham et al.l ( l2009h observed a silicate feature with a crys- 
talline peak, in contrast to the pre-outburst silicate feature which 
was similar to the interstellar medium spectrum dominated by 
amorphous grains. They concluded that thermal annealing in the 
surface layer of the proto-planetary disk was the mechanism for 
this change in spectral features. Unfortunately, we have no pre- 
outburst spectrum of V 1 1 1 8 Ori before its outburst to compare 
with the outburst spectra, only a spectrum taken 2 years after the 
end of the outburst and 3.5 years after our last outburst spec- 
trum in March 2005. If the pre-outburst spectrum was indeed 
typical of amorphous silicates, it means that the time scale for 
disappearance of the crystalline features must be longer than 2 
years. Abraham et alj ( 1200 9) proposed that similar crystalliza- 
tion must have occurred during the previous 1955-1956 outburst 
of EX Lup, and since the 2005 quiescent spectrum did not show 
crystalline features, the removal time scale must have been less 
than 50 years. On the other hand, in VI 118 Ori, it is also pos- 
sible that the unchanging silicate feature in the three spectra is 
typical of the quiescent silicate feature in the young star, per- 
haps an indication that the disk is more evolved. In any case, 
M-type young stars typic ally show flat silicate features (e.g., 
iKessler-Silacci et al.ll2005[ 12006). Finally, we note that the IRS 
data do not show strong evidence of the forsterite crystalline fea- 



ture at 16 ^m, in contrast with the feature observed in EX Lup 
([Abraham et al.l2009h . This is consistent with the lack of a sharp 
peak at 10 fj,m and of a shoulder at 11.3 /im. 

5.2. Color indices 

We have calculated color indices from our optical and near- 
infrared data, and looked into their evolution during the out- 
burst. Time intervals were defined in the pre-outburst epoch 
(1; JD < 2453200), initial phase of rise (2; 2453200 <JD < 
2453390), second phase of rise (3; 2453390 <JD < 2453450), 
plateau (4; 2453450 <JD < 2453650), decay (5; 2453650 <JD 
< 2453900), and post-outburst (6; 2453900 <JD < 2454200). 

Figure [10] shows the color-magnitude diagrams, whereas 
Figure [TT] shows the color-color diagrams, together with an 
Ay = 1 mag re ddening vector using the reddening law of 
iRieke et all ([1985'), the loci of the unreddened main sequence 
and giant stars (Bessell & Brett 1988), and the loci of classical 
T Tauri stars (Mever et alj|1997l) . Compared to the pre-outburst 
color indices of VI 1 18 Ori, the optical color indices are bluer in 
the early phases of the outburst (except for B — V which shows a 
peculiar pattern, probably because the B magnitude is very sen- 
sitive to the evolution of hotspot emission), and gradually return 
to the pre-outburst values in the decaying phase of the outburst. 
In the near infrared, an apparent blueing occurs as well in the ini- 
tial phases of the outburst with a return to "normal" colors in the 
decay phase. The above results are consistent with the hypoth- 
esis of Lorenzetti et al. (2007) that blueing must have occurred 
due to a hotter temperature component during the initial phase 
of the outburst (their data were taken only in the decaying phase 
and showed reddening). The variations in the color-color dia- 
grams also show no trend related to extinction variations, as sug- 
gested by Lorenzetti et al. (2007) and Lorenzetti et al. (2009). 
There is, however, a larger scatter in the infrared color-color dia- 
grams, especially in J—H versus H — K. Since the stellar pho- 
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Table 4. Spectral fits and quantile properties. 



Parameter Sep 2002 

Satellite Chandra 
Spectral analysis: 

Net counts" 150.9 

Scaled background counts" 27.1 

Exposure (ks)" 47.85 

A^H (lO^i cm-2) 2.7lJ;^ 

T{MK) 25.lt|g 

EM(10^^cm-3) lOloi 

Z/Zq :=0.17 

Fx (10-^"* ergscm-2 s-i) 3.4l°-^ 

Lx (10=^0 ergs s"!) 1.2+^;J 

Quantile analysis: 

Energy range 0.2 — 6.0 

Net counts'^ 150.44 

Scaled background counts'^ 26.56 

Q25% 0.17 ±0.01 

Q50% 0.22 ±0.01 

Q75% 0.31 ± 0.03 

Ego%(keV) 1.45 ±0.06 

X = log[Q50%/(l - Q50%)] ■ ■ ■ ■ -0.56t°-E;^ 

y = 3 X Q25%/Q75% l-63iE;-i^ 

T^B (MK) Ri 20 

NYi.xy (lO^i cm-2) ~ 3 



Jan 2005 



Feb 2005 



Mar 2005 



Sep 2005 



Chandra 

10.8 
0.2 
4.66 

1 4+3-6 



46 
0.4 

2.0 
0.5 



28 
+0.4 
0.1 

0.17 

1.5 
1.1 
+0.4 
-0.3 



0.2 - 4.0 
10.79 
0.21 
0.24 ±0.07 
0.33 ±0.07 
0.42 ±0.14 

1.44 ±0.25 

1+0.12 
-0.14 



-0.31"' 



1 TO+O-Si 

^- '"-0.61 
« 12 



XMM-Newton 

285.5/105.4/103.2 
224.5/56.6/52.8 
18.34/21.87/21.91 



4.3 
7.7 
2.9 
0.17 
2.9 
2.5 



+ 1.3 

1.1 

+ 1.3 
0.8 
+ 1.0 
1.0 
+0.20 
0.06 
+0.2 
0.2 
+0.2 
0.2 



0.2 - 4.0 
286.45 
222.55 
0.15 ±0.01 
0.20 ±0.01 
0.29 ±0.01 
0.95 ± 0.04 
-0 61+° °^ 

'^■°^-0.03 

^ 12 
Si 1 



XMM-Newton 

{47.1)/190. 4/108.0 
(698. 9)/325. 6/331.0 
{7.83)/19.96/19.22 

1.4^ ' 
15.9^ 
2.1-' 



1+0.8 
-0.5 
3+3.9 
-2.0 
1+0.2 
-0.2 
: 0.17 
3+0.7 
-0.7 
9 1+0.2 
^•-^-0.2 



5.9"* 



0.2 - 4.0 
196.34 
425.66 
0.15 ± 0.02 
0.25 ± 0.02 
0.34 ± 0.03 
1.14 ±0.07 
-0.48_o 04 

1 ,n+0.19 

a 20 



XMM-Newton 

(12.4) 
(227.6) 
(13.03) 



1.2^ 



,+1.2b 
-1.0 
0.04 - 0.12'' 



0.2 - 4.0 
9.75 
228.25 
0.0 ±0.5 
0.3 ± 1.5 
0.9 ± 1.3 
1.2 ±4.0 



-0.4 
0.05^ 



+ 00 



; O.i 



Jan 2006 Feb 2006 Mar 2006 Apr 2006 Jul 2006 Dec 2007 

Satellite Chandra Chandra XMM-Newton Chandra Chandra Chandra 
Spectral analysis: 

Net counts" 30.5 10.5 623.4/183.8/232.4 36.4 36.3 37.9 

Scaled background counts" 0.5 0.5 998.6/215.2/209.6 0.6 0.7 1.1 

Exposure (ks)" 27.82 27.11 78.38/90.90/90.95 29.57 29.57 34.73 

iVn (10^1 cm-2) 3^+2.4 3 5+4.3 ^ gl?;? 4 9+5.2 7 3+2.2 ^^^+1.6 

T(MK) 46.81^0 8.7l^-^ H-Sljg / 89^43 15. 9t^-^ 6.5^2? 19.9t^-* 

EM(10^3^m-3) 0.18l° i2 O.W+J;\l / 0.51+^;^ 0.35t° «? 0.20l° 08 

Z/Zq :=0.17 :=0.17 ':= 0.17 ' := 0.17 := 0.17 := 0.17 

Fx (10-1* ergs cm-2 s-i) o.9lt« i« 0.13l» °i ^i.l+Jil 0.63l° lJ 0.47« «? ^-^^^-IZ 

Lx (10^0 ergs s-i) 0.271^^ O.iol^^ 1.16«;«^ 0.361^^ 0.9«;? 0.201^^ 

Quantile analysis: 

Energy range 0.2-4.0 0.2-4.0 0.3- 10.0 0.2-4.0 0.2-4.0 0.2 -4.0 

Netcounts<= 30.48 10.53 626.9 36.37 36.32 37.91 

Scaled background counts'^ 0.52 0.47 985.1 0.63 0.68 1.09 

Q25% 0.28 ±0.04 0.18 ±0.03 0.081 ± 0.010 0.23 ± 0.02 0.18 ± 0.01 0.22 ±0.03 

Q50% 0.35 ±0.04 0.24 ±0.09 0.145 ±0.006 0.31 ± 0.03 0.24 ± 0.04 0.28 ± 0.02 

Q75% 0.48 ±0.08 0.46 ±0.19 0.255 ± 0.019 0.40 ± 0.06 0.34 ± 0.03 0.40 ± 0.06 

F5o% (keV) 1.52 ±0.15 1.12 ±0.35 1.709 ± 0.06 1.39 ±0.09 1.11 ± 0.14 1.28 ±0.09 

x = log[Q50%/(l-Q50%)]-.-- -0.27to«^ -0.50lg-^o -0.77tn2 -0.34l°-°5 -0-50«-«« -OAQ+_°oll 

y = 3 X Q25%/Q75% l-74lE!fg 1.19l2f2 0-95ton l-esl^lz ^■^^tl 'S, l-GSlSfg 

T^a (MK) Si 13 35 R:i 90 Ri 10 ^15 Ri 16 

NTi,xy (10^^ cm-2) -10 ^ 0.2 2.5 ^ 9.0 1.5 ^ 3.5 

The uncertainties are based on 68% Bayesian confidence ranges. X-ray luminosity and absorbed X-ray flux in the 0.1 — 10 keV range, assuming d = 400 pc. 
^ ¥or XMM-Newton observations, values refer to the EPIC pn, MOSl, and MOS2, respectively (pn only for September 2005). The values in parentheses refer to 

detectors not used in the spectral fits. 
^ Estimates based on models for February 2005 and March 2006. 

Based on pn data in general, except for March 2005 for which MOSl data only were used (similar values and results were obtained for MOS2 data). 



tospheric emission does not vary during the outburst, we need to 
understand the variations in color indices and magnitudes dur- 
ing the outburst in the context of a star + disk + hot spot model. 
iLorenzetti et alj ( l2007b suggested that the polarization found in 
the / band of the V 1 1 1 8 Ori outburst was intrinsic and that it was 
probably due to the spotted and magnetized stellar photosphere. 
As we will show below, the photometric changes can indeed be 
explained by changes in the fluxes of the disk thermal emission 



and from a hot spot as the amount of mass falling from the disk 
onto the stellar photosphere increases during the outburst. We 
note also that during the outburst, any increase in hotspot emis- 
sion will further increase the irradiation of the disk. The follow- 
ing section aims to understand the multi-wavelength photometry 
and spectroscopy observed during the outburst of VI 118 Ori by 
comparing the data with young star disk models. 
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Fig. 8. Color-color diagrams based on the (550% (= m), Q25%, and Q75% quantiles. The grid is in plasma temperature (kT in keV) 
and hydrogen column density (N^ in 10^^ cm^^). The serendipitous, pre-outburst observation of September 2002 is shown at the 
top, the post-outburst version of December 2007 is shown at the bottom. The 9 outburst observations in between are ordered in time 
from left to right and top to bottom, with January 2005 at the top left and July 2006 at the bottom right. 



5.3. Spectral energy distribution 

Figure [12] shows the observed SED for VI 118 Ori. We have 
shown as small colored dots the data for the different phases 
(as defined at the beginning of section 15.21 and in the caption 



of Fig. [TOb . and the pre- and post-outburst photometry as red 
crosses and black stars, respectively. The IRS low-resolution 
spectra are also shown for the outburst phases (violet) and post- 
outburst phase (black). For clarity, we do not show the high- 
resolution spectra. To emphasize the expected stellar photo- 




Fig. 9. Spitzer IRS SL spectra of VI 118 Ori for the three epochs (left panel). The dotted line shows the fit to the underlying 
continuum in the [6-8,13-14.5] fim range. The right panels show the normalized flux at the three epochs. 






H-K J-K V-K 

Fig. 10. Magnitude versus color indices for the different time intervals (red: #1, lime: #2, violet: #3, blue: #4, green: #5, black: #6). 



spheric contribution, we have added a stellar atmosphere model 
from the PHOENIX databas^ (e.g., Hauschi ldt e t al. 1999) for 
a stai- with T^s = 3600 K, log 5 = 4.0, and Z = 
(STARdusty2000 models ), and reddened it with Ay = 1.5 mag 
and the reddening law of lFitzpatrickl (Il999t) . Indeed, the hydro- 



ally staying in the range 1 ~ 5 x 10^^ cm ^, which corresponds 
to Ay in the range Ay = l.Tl^pgy mag for Ry = 3.1 (Galactic 
value) or Ay 
ment). 



1.44lQ4g mag for Ry — 5.5 (dusty environ- 



gen column density derived from X-rays is 2.7^0 g 
in 2002, and does not change much during the outburst, gener- 



We calculated the median magnitude differences from 0.44 
to 8.0 /im (we have no U photometry during the outburst), rela- 
tive to the median magnitude in quiescence (we used both pre- 
|http : / /www ■ hs . uni-hamburg . de/EN/For /ThA/phoeni;Qtahifek aadnT3i) St-outburst magnitude; Fig. [T3]l. There is cleai" 
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Fi g. 11. C olor index diagrams for the same time intervals as in Fig.[TO] An Ay = 1 mag vector is added using the reddening law 
of lRieke e t al. ( 1985). The loci of unreddened main sequence and giant stars ( iBessell & Brett 1988) are shown for the near-infrared 
panels as solid and dotted curves, respectively. Finally the loci of classical T Tauri stars jMever et al.ll 19971) are shown as a dashed 
line in the (J — H) vs. (H — K) diagram. 



evidence that the optical emission brightened much more than 
the infrared emission, with a peak in the V band. We also ob- 
serve that the BVRI decayed faster than the JHK bands, prob- 
ably indicating a mixture of different emission components in 
the near-infrared bands (e.g., disk and hot spot). At long wave- 
lengths (A > 2 iim), the increase in flux was similar at all wave- 
lengths. 

To quantify better the star-ndisk properties of VI 1 18 Ori be- 
fore and during its outburst, w e have then used the code de- 
scribed in lWhitnev et alj (l2003a"b). version 2009022^ We have 
run the code by fixing certain parameters and leaving others 
free to vary. First, the models did not include an envelope, but 
only an ambient medium with density of 10^^^ g cm^^. The 
contribution of this medium is generally negligible. The exclu- 
sion of the envelope is motivated by the fact that it would con- 
tribute mostly at long wavelengths where we have no data. The 
IRS data also do not show strong evidence of an envelope (that 
could not be attributed to improper background subtraction). The 
fixed parameters were the stellar radius (i?^ ~ 1.29i?0), mass 
(M^ = O.4IM0), and effective temperature (Tcff = 3600 K; see 
ij2l we used as input a lower spectral resolution version of the 
same PHOENIX stellar atmosphere model as above), the disk 



^ |http : / / gemelli ■ Colorado ■ edu/ -bwhitney/ codes/| 



mass (5Afjup) and outer radius (100 AU), the disk flaring power 
(cx Rl^, (3 = 1.1) and the disk scale height (O.Oli?,^ at i?*; the 
height is thus 2.6 AU at 100 AU), and the disk density profile 
(oc i?^", a = 2.1). For the dust opacity, we used a model of 
grains as large as 1 mm for the denser regions of the disk (see 
Wood et al. 2002 for HH 30), but used an interstellar medium 
grain model for the less dense regions of the disk (i.e., below a 
threshold set at 3.34 x 10^^'* cm^'^). The disk mass and outer 
radius were not constrained from our optical and infrared data0. 

Changes in their values would also not affect much the mod- 
els at optical-infrared wavelengths. We have left the mass accre- 
tion rate and inner disk radius (in units of the dust destruction 
radius, i?sub) free to vary. The model also includes the possi- 
bility for material to fall from the disk onto a hot spot by fol- 
lowing the magnetic field lines. We have kept the truncation ra- 
dius to 5 AU, but left the fractional area, fs (ratio of the shock 



We have obtained data with the IRAM 30m MAMB02 bolometer at 
1.2 mm in late November 2009. VI 1 18 Ori is detected at 1.8± 0.4 mjy. 
Assuming an optical thin regime, an opacity of k„ = 0.02 cm^ g~^, us- 
ing a gas-to-dust ratio of 100, we obtain a disk mass of 2.45 ± 0.5Af jup 
for a characteristic temperature of 20 K. The disk mass increases to 
6.9 ± l.SAfjup for Tc = 10 K. Although the determination of the exact 
disk mass is out of scope in this paper, the derived estimates indicate 
that the assumed disk mass in our model is realistic. 
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Fig. 12. Compilation of the optical and infrared photometry taken before, during, and just after the outburst (see text and caption of 
Fig. [To] for details about colors). The IRS SL and LL spectra and a reddened (Ay = 1.5 mag) stellar atmosphere model are also 
shown. 



area to the star area; see ICalvet&Gullbring|ll998h . free to vary. 
Finally, the code samples the model at 10 different inclination 
angles from cosi = 0.05 to 0.95 every 0.10, i.e., from 87° 
to 18°. Thus, we have four different parameters that can vary 

([C0SZ,M, /s,i?min])- 

Note that the choice of a and f3 is typical of those observed in 
T Tauri stars, and correspond to values that seem to fit the VI 1 18 
Ori data. The power law indices are related in the sense that if 
the disk surface density is E oc RP, then p = f3 — a. A value 
of /3 = 1.1 is close to the theoretical values for irradiated and a 
disks (/? = 9/8 = 1.125), and, with a = 2.1, p = /3-a = -1.0, 
i.e., typical for a steady optically thick disk. Note that a fit to 
the IRAC and MIPS fluxes in quiescence gives \F\ cx A^°^^, 
i.e., consistent with the commonly observed XFx oc X~^^^ for 
irradiated disks. Assuming that the temperature follows a power 
law T oc i?^"?, and if XF\ oc A*, both indices are related by 
s = 2/(7 - 4. Thus, we have q = 0.6, i.e., T oc R-^/'^. 

There is no pre-outburst IRS spectrum available to determine 
the shape of the silicate feature before the outburst (emission or 
absorption). However, we have used the post-outburst IRS spec- 
tra, scaled by a factor of 0.7 to match the IRAC photometry (see 
Fig.|9j. For the photometry, we have used the median of the pre- 
outburst (interval #1) and post-outburst (interval #6) photometry. 
For the outburst data, we have used the plateau (#4) peak pho- 
tometric data with the second epoch outburst IRS SL spectrum 



to characterize the SED at "peak". The IRS spectrum akeady 
matches well the 8 /um flux without correction and shows the 
silicate feature in emission. 

We ran the code with 5,000,000 "photons" or energy pack- 
ets for each model. This results in SED models with reasonably 
good signal-to-noise ratios over the wavelength range covered by 
our data. The models were then reddened with Ay ~ 15 mag 
(assuming ISM opacities). For each reaUzation we have calcu- 
lated a reduced Xp value 



N 



N ^ 

1 



A(logioF,(A,)) 



(1) 



where F^{Xi) are the data flux densities, 
Mi,{Xi,p) are the model flux densities (for pa- 
rameter vector p = [cos i, M, fs, Rmin]) in- 
terpolated at the wavelength of the data points 
(0.44, 0.55, 0.71, 0.97, 1.25, 1.60, 2.22, 3.6, 4.5, 5.8, 8.0, 24.0 
for the photometry, and 5.5, 6.0, 6.5, . . . , 14.0 /im for the IRS 
- we only used the SL data). For the photometric flux density 
uncertainties, we have used a fixed A(logj^Q i^i,(Ai)) = 0.05, 
i.e., about 10% relative uncertainties for all data points (a value 
larger than many actual uncertainties on the measured values) 
to account for flux variations observed before and at the peak of 
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Wavelength (|im) 

Fig. 13. Median magnitude difference relative to quiescence as a 
function of wavelength. The same time intervals as in Fig.fTOlare 
used for the color coding (lime filled circles: #2, violet triangles: 
#3, blue diamonds: #4, green stars: #5). In the initial phase of 
the outburst, a strong optical component developed, with a peak 
in the V band. 
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Fig. 14. Probability distribution of the parameters 
([cos i, M, fs, Rmin]) for the quiescent (solid) and outburst 
(dotted) cases. 



the outburst. We used A(log^o Fvi^)) = 0.025 instead for the 
IRS SL data to put more weight on the silicate feature. 

Figure [14] shows the normalized model probabilities (calcu- 

lated as pp = e ^p) for the quiescent case (solid), and for the 
outburst case (dotted). The distributions of inclination angles are 
slightly different for the quiescent and peak of outburst cases, but 
the distributions overlap in the range cosi = 0.1 — 0.3, with a 
maximum probability (obtained by multiplying the probability 
distributions) at cosi = 0.25, i.e., i « 76°. The mass accretion 



rate distributions are sharply peaked: the most probable values 
of M are lO"*^'^ Mq yr^^ and 10"'^ ° M© yr^^ in quiescence 
and and at the peak of the outburst, respectively, indicating an 
increase by a factor of 4. The coverage factor seems to be small 
(fs < 0.0001) outside the outburst while the simulations indi- 
cate high values (0.1 < fs < 1) at the peak of the outburst. Note 
that this parameter is essentially constrained by the optical and 
near-infrared photometry. The minimum disk radius has likely 
values of i?inin = (10° * - 10°-^)i?sub, with the most probable 
value at i?niin ~ 4i?sub, whereas the values at the peak of the 
outburst are lower and most likely equal to the sublimation ra- 
dius. The best models give sublimation radii of about 0.10 AU 
and 0.20 AU in quiescence and at the peak of the outburst, im- 
plying inner disk radii of 0.25 — 0.6 AU (most probable value 
of 0.4 AU) and 0.2 AU, respectively. This result suggests that 
the inner disk radius (determined from the infrared data arising 
from dust emission) probably diminished during the outburst, 
although we emphasize that quiescent values of 0.2 AU are pos- 
sible, though less probable from our SED modeling. 

Figure [15] shows the best fits to the "quiescent" SED and to 
the peak of outburst SED, respectively. The total (bolometric) 
luminosity (including all components) varied from about 2.0Lq 
in quiescence (with a contribution of the stellar luminosity, L — 
AnR^aT^ff, of O.25L0 and of the thermal disk luminosity of 
0.9^0) to « 7.4^0 at the peak of the outbursQ (the thermal 
disk luminosity increasing to about I.7L0). 

Thus the bolometric luminosity is dominated by the stellar 
direct emission, i.e., the stellar photospheric emission and the 
hotspot emission (the latter in particular during the outburst). We 
note, however, that we would need data at shorter wavelengths 
than B to constrain better the hotspot emission during the out- 
burst. 



6. Discussion 

Our multi-wavelength campaign allows us to obtain several 
pieces of information on the outburst of VI 118 Ori, as seen in 
the optical and infrared, and in the X-rays. Overall, the evolu- 
tion of the X-ray flux during the outburst followed the evolu- 
tion of the optical and near-infrared fluxes. The sampling of the 
X-ray observation does not allow us to determine whether the 
increase in X-ray flux in March 2006 is related to the short- 
term increase in optical/infrared flux observed about 45 days 
earlier. The thermal properties of V 1 1 1 8 Ori showed evidence 
of a plasma with temperature of a few MK to a few tens of MK. 
In our initial paper, we reported a change in the spectral proper- 
ties from a predominantly hot (20 MK) plasma before the out- 
burst (and possibly at the ver y beginning) to a lo wer temperature 
(8 MK) in February 2005 (Auda rd et all 2005b). Our reanalysis 
confirms this finding and a gradual return to larger values in the 
later phases of the outburst. However, the low signal-to-noise ra- 
tio of the X-ray spectra taken during our monitoring campaign 
makes it difficult to constrain accurately the plasma temperature 
(in contrast, the February and March 2005 observations had a 
much better signal-to-noise ratio). The alternative quantile anal- 



' iLorenzetti et al.l i2006l) quote a peak luminosity of 25.4L0; how- 
ever, they used the peak magnitudes of a previous outburst to build the 
outburst peak SED and to derive the luminosity. The previous outburst 
had a peak magnitude of J = 8.5, in contrast to the 2005-2006 outburst 
that peaked at J = 10.6. In addition, our method derives the lumi- 
nosity by ma tching the observ ed SED with a detailed radiative transfer 
model, while iLorenzetti et al.l |2006) integrated the flux densities. The 
two above issues may explain the difference in derived luminosities. 
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Fig. 15. The best-fit models to the "quiescent" SED (left) and to the peak outburst SED (right). The direct stellar contribution 
(including emission due to the hot spot) is shown as a dotted line (the Ay = 1.5 mag reddened stellar atmosphere model is shown 
as a dash-dotted line), the thermal disk spectrum as a dashed line. 



ysis showed evidence of less variations in the plasma tempera- 
ture. Nevertheless, it is consistent with the detected change in 
temperature. While V1647 Ori showed a strong increas e in X- 
ray fl ux during the outburst (Kastner et al. 2004; Gross o et alj 
12005b . VI 118 Ori only showed moderate X -ray flux increase, 
and possibly an eventual decrease in flux after the outburst, if 
we assume that the 2002 flux level was representative of the pre- 
outburst X-ray flux. In this case, the impact of the mass infall 
onto the stellar magnetosphere must have significantly affected 
the latter, leaving it in a lower X-ray flux state. It is, however, 
also possible that the 2002 flux was observed by chance in a 
high quiescent state, and that the 2007 X-ray flux is closer to 
the pre-outburst flux level. In any case, the changes in the X-ray 
fluxes during the outburst, and possibly the temperature changes 
all indicate a strong interplay between the material falling from 
the disk due to the increase in mass accretion rate and the stellar 
magnetosphere. 

From the optical an d infrared SEP data a nd Monte-Carlo 
models using the code of IWhitnev et al.l (|2003a''b'), there is clear 
evidence that the thermal disk emission dominates in the mid- 
infrared in quiescence and during the outburst, while emission 
from a hotspot and reprocessed emission dominates the outburst 
SED in the optical and near-infrared, and contributes to some ex- 
tent before the outburst. However, we note that this conclusion 
depends on our choice of using disk models that stop at the sub- 
limation radius, i.e., there is no disk contribution at shorter radii, 
where the gas is optically thin. This could contribute to the near- 
infrared fluxes, but it should not contribute significantly to the 
optical fluxes. Therefore, we believe that the addition of hotspot 
emission is necessary to fit the blue side of the VI 1 18 Ori SED, 
during the outburst and probably before the outburst. 

The increase in mass accretion rate during the outburst leads 
to the change in inner disk radius which moves closer to the 
star's surface. Such an effect could explain the changes in the 
X-ray emission (see above). From our simulations, we have de- 
termined that the mass accretion rate increased during the out- 
burst: consequently, the dust sublimation radius increased from 
« 17 i?* = 0.1 AU to w 32 i?,^ = 0.2 AU, whereas the inner 
disk radius (based on the infrared dust emission) moved from 
0.4 AU in quiescence to 0.2 AU in outburst (although the op- 
tically thin gas disk could be closer to the star), and the inner 
disk radius is consistent with the dust sublimation radius at the 
peak of the outburst. We emphasize that such estimates of the 



inner disk radius are biased due to our use of dust thermal emis- 
sion for the disk. The optically thin gas disk is difficult to de- 
tect and we would require gas tracers, e.g., H2 or CO bandheads 
to properly determine the inner extent of the disk. In any case, 
the strong increase in fractional area of the hotspot suggests that 
large amounts of matter fell from the disk onto the star during the 
outburst, incidentally interacting with the stellar magnetosphere. 

Using their near-infrared spectrum taken in September 2005 
(i.e., a bout at the end of the plateau phase) and a wind loss 
model. iLorenzetti et al.l ( l2006h derived a mass loss rat^l of about 
4 X 10~^ Mq yr~^ from Hi. Furthermore, little or no ab- 
sorption was reported (Ay < 2), with no evidence of varia- 
tions during the outburst. They also derived a mass loss rate 
of (3 — 8) X 10~^ Mq yr~^ from the CO overtone emis- 
sion, indicating an ioni zation fraction of 0.1 — 0.2. Furthermore, 
ILorenzetti et al.l ([2007^ observed VI 1 18 Ori in September 2006, 
i.e., well after the outburst, and detected no emission line, while 
the 2.3 /im CO band was seen in absorption. Evidence of wind 
emissi on in Ha (P Cyg profile) was also reported by iHerbid 
(I2008h in November 2005, in contrast with the more symmet- 
ric profile observed after the outburst in December 2006. It is, 
therefore, likely that high wind loss rates occurred transiently 
due to the instability in the disk. We note that disk emission may 
provide a better fit to C O bands in young sta rs than mass loss 
emission in general (e.g.. lChandler et aLlll995h . but this may not 
be the case during the outburst of a young star such as VI 118 
Ori: our SED models including emission due to matter falling 
from the accretion disk onto the stellar photosphere clearly show 
that the outburst SED at wavelengths below 3/im is dominated 
by the hotspot emission, not by disk thermal emission (while it 
remains negligible outside the outburst and stellar photospheric 
emission dominates). Wind loss emission may thus be present in 
the optical and near-infrared. 

We can compare the mass loss rates derived for the neu- 
tral gas by Lorenzetti et al. (2006) to our derived mass accretion 
rates: the SED modeling suggests a variation in the mass accre- 
tion rate from about 2.5x10"'' Mq yr^i to 1.0 xlO"'' Mq yr~\ 
i.e., these mass accretion rates are comparable to the mass loss 



^ ILorenzetti et al.l l f2006l) derive the mass loss rate from the wind 
model of ICarj ("1989) and his Figure 8. While their derived mass loss 
rate is accurate, they incorrectly give a CO luminosity of 6.3 x 10~^ Lq, 
while it should have been 7.9 x lO~*I/0 based on the quoted line flux 
(D. Lorenzetti 2009, priv. comm.). 
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rates determined by iLorenzetti et alj (l2006h . This is somewhat 
surprising, since mass outflow rates are expected to be about 0.1 
times lower than mass accretion r ates ( Hartigan et al. 19 951 with 
revise d mass accretion rates, see iGuUbring et al.lll998l:IEdwardsl 
120091) . Our mass accretion rate in quiescence is larger by about 
a factor of 10 compared t o those derived for cla ssical T Tau stars 
(« IQ-^Mq yr-i; e.g.. lGullbring et al.l ll998'). and lower by a 
factor of 10 than the r ates derived for emb edded Class 1 stars 
(« 5 X IQ-^Mo vr-^: lKenvon et al.lll993l) . Therefore, our de- 
rived mass accretion rates are consistent with placing VI 1 18 Ori 
in a category between Class 1 and CTTS stars. We note also that 
the IRAC color indices of VI 118 Ori in quiescence are [3.6]- 
[4.51=0.65-0.75 and [5.8]-[8.0]=0.65-0.7 which places it for- 
mally in the region of Class II sources, b ut also close to the re- 
gion for Class O/I (e.g.. lAllen et alj|2004 . Furthermore, the K- 
[3.6] color of 1.52 in "quiescence" also place the star as a strong 
accretor. In contrast, U — V = —0.05 which formally would 
(incorrectly) place VI 118 Ori as a w eak-lined T Tauri star In 
fact, using the lGullbring et al.l (Il998h method to derive mass ac- 
cretion rates from the t/-band luminosity excess, we find a mass 
accretion rate of 4.1 x 10~^ AIq yr^^, much smaller than val- 
ues derived by SED modeling, and also inconsistent with the 
strong infrared excess. While our derived mass accretion rates 
are model-dependent, we are confident that our derived mass ac- 
cretion rates are realistic, despite the difficulty to derive mass 
accretion rates from infrared excess (e.g., Kenyon & Hartmann 
Il987h . 

Overall, the data provide ample evidence that a significant 
portion of matter fell from the disk onto the young star due to 
the increase in mass accretion rate. Over a duration of about 
400 days, we estimate about 5 x 10^'' Mq — O.ISA/q was de- 
posited (using an average of half the peak mass accretion rate), 
i.e., about twice the amount that would have been deposited over 
the same time range if using the quiescent mass accretion rate 
The hotspot covering factor we derive at the peak of the outburst 
is also relatively large (fs > 0.1), suggesting that a large frac- 
tion of the stellar magnetosphere may also have been influenced 
by the increase in mass accretion rate, which in turn supports the 
hypothesis that the stellar magnetosphere was modified during 
the outb urst and the X-ray pr operties of the coronal plasma have 
changed lAudard et alj (l2005b.) . Our SED modeling of the dust 
thermal emission suggests a change of the inner dust disk radius 
from about 0.4 AU in quiescence to 0.2 AU at the peak of the 
outburst (the optically thin gas could reach down to 0.2 AU or 
below before the outburst, although it would be unclear why dust 
would not go down to the sublimation radius, unless, e.g., disk 
clearing by a giant planet or significant grain growth making the 
dust optically thin would be present in the inner portion of the 
disk). Together with the increase in fs, this provides evidence 
that the disk closed in during the outburst, disrupted the mag- 
netosphere, and matter fell onto the stellar photosphere, produc- 
ing strong hotspot emission that dominated the optical and near- 
infrared in outburst (while the thermal disk emission increased 
by a factor of 4 and dominated the mid-infrared emission). 

While color-color and color-magnitude diagrams can be used 
in young stars to determine the level of reddening, such diagrams 
might be difficult to use to measure any reddening during the 
outburst of a young star: in addition to reddening, the whole SED 
of a young star changes dramatically, increasing in flux, but also 
in shape: in the case of VI 1 18 Ori, the optical and near-infrared 
fluxes increased much more than the mid-infrared spectrum. The 
former are dominated by stellar and hotspot emission, while the 
latter is dominated by thermal disk emission. Color-color dia- 
grams show a clear blueing of the colors (Fig. fTTT i. not due to 



clearing of any enshrouding material, but due to the strong in- 
crease in flux below 1.5/im (Fig. [T3). Such multi-wavelength 
data show the importance that a good coverage of the SED is 
necessary to understand the disk-star interactions in young stars. 
Furthermore, Figures[TO]and[TT|indicate that VI 1 18 Ori returned 
(within uncertainties) to its "initial", pre-outburst conditions. 



7. Conclusions 

We have followed the young accreting star VI 118 Ori before, 
during, and after its outburst in 2005-2006 with photometry in 
U BVRIJHK and in X-rays with XMM-Newton and Chandra. 
We have further obtained Spitzer photometry and spectroscopy 
during the outburst and after the outburst. 

The optical and infrared data showed significant variations: 
the data below 1 jjm brightened by as much as about 3 — 4 mag, 
while the infrared data brightened only up to « 2 mag. Monte- 
Carlo simulations of a star-ndisk-nhotspot model suggested that 
the optical data were dominated by hotspot emission at the peak 
of the outburst (while the stellar photosphere dominated before 
the outburst), while thermal disk emission from optically thick 
dust dominated the infrared. The SED analysis showed an in- 
crease in mass accretion rate from M — 2.5 x 10~^ yr~^ 
in quiescence to M = 1.0 x 10~^ Mq yr^^ at the peak of the 
outburst, and an increase of the hotspot fractional coverage area 
from Js < 10~^ to fs > 0.1. The SED modeling suggests that 
the dust sublimation radius increased from 0.1 to 0.2 AU. Based 
on the thermal dust emission, the inner disk radius moved from 
about 4 times the dust sublimation radius (0.4 AU) to 0.2 AU 
{— i?sub) at the peak of the outburst. But the gas inner disk radius 
could extend below the dust sublimation radius. The inclination 
angle of the disk is constrained to values of cosi — 0.1 — 0.3 
with a peak at 0.25, i.e., z « 76°. 

The Spitzer IRS spectra showed the silicate feature at 10 fim 
in emission during and after the outburst, with no apparent vari- 
ation in shape or flux, indicating that the outburst had little im- 
pact on the optically thin region where the silicate emission is 
originating. Emission lines were detected in the high-resolution 
module data although the background subtraction proved to be 
difficult, especially for the outburst spectra. Follow-up observa- 
tions, e.g., from the ground should confirm the detection of lines 
of particular interest, e.g., [Ne ll], [S III], and H2 - 0. 

The initial X-ray observations (Jan-Mar 05) indicated a 
change in the thermal structure of the coron al spectrum of V 1 1 1 8 
Ori with httle change in the X-ray flux ( Audard et al.ll2005bl) . 
while fading in the X-rays w as already detected in September 

2005 (Lorenz etti et al.l 12006). Our continuing monitoring with 
XMM-Newton and Chandra provided the complete picture, con- 
firmed the low plasma temperature measured in February 2005, 
and showed that VI 1 18 Ori's X-ray flux was strongly correlated 
with the optical and infrared fluxes, although intrinsic variabil- 
ity (due to magnetic activity) was also observed. Compared to 
the 2002 flux, the increase in X-ray flux was not large during 
the outburst and the post-outburst X-ray observation in late 2007 
showed a lower flux level than observed serendipitously in 2002, 
which could either indicate that the outburst significantly im- 
pacted the stellar magnetosphere to leave it in a lower state than 
before the outburst, or that the 2002 observation caught VI 118 
Ori in a quiescent high state. The X-ray spectra showed evidence 
of thermal variations of the corona during the outburst, but the 
signal-to-noise ratio of our data was not large (except in March 

2006 during which we could detect a second temperature com- 
ponent, probably due to the increase in flux at the end of the 
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observation, likely contaminated by a flare). After the apparent 
change to a cooler plasma temperature in the early phase of the 
outburst, we argue that the thermal properties of VI 118 Ori's 
corona returned to its "normal", hotter values before the end of 
the outburst, but that the continuing mass infall affected signifi- 
cantly the magnetospheric configuration. 

Our observations have shown the interplay between disk and 
stellar magnetosphere in young accreting stars, as demonstrated 
by the changes in the X-ray plasma emission due to the strong 
increase in mass infall from the disk onto the stellar photosphere. 
While our data do not directly help to determine the origin of the 
outbursts (thermal disk instabilities or instabilities due to close 
companions or giant planets in the disk), it is interesting that 
iReipurth etaTI ( l2007h identified V1118 Ori as a close visual bi- 
nary. Future observations should, ideally, take into account this 
binarity to determine which star is actually outbursting. 
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Appendix A: Details on tlie SMARTS data 
reduction 

The ANDICAM optical detector is a Fairchild 447 2048 x 2048 
CCD with 30/im pixels. It is used in 2 x2 binning mode, yielding 
a nominal plate scale of 0.369 arcsec pixel^^ and about a 6' x 6' 
field. Readnoise is 6.5 e^ RMS; gain is 2.3 e^/DN. It is read 
with a single amplifier. Filters are Johnson B and V and Kron- 
Cousins R and /. We generally took 3 exposures in B and V, 
and 2 images in Rc and Ic, and co-added these prior to analysis. 
The IR channel detector is a Rockwell 1024^ HgCdTe "Hawaii" 
array with 18 /im pixels binned 2x2 to yield a 0'.'274 pixel^^ 
plate scale. The IR channel field is 2'A x 2'A. The filters we used 
are standard CIT/CTIO JHK. We flatten the images using dome 
flats obtained every 3 nights. We observe using a 3-point dither 
pattern and a 40" throw, and use the median of the unshifted 
images as the local sky. 

A.1. Differential Photometry 

We initially determine relative magnitudes using differential 
photometry. In the optical, we selected the 1 1 brightest stars in 
the « 6' field of view (Table lA.ll V372 Ori is overexposed). 
Instrumental magnitudes are determined for the target and all 
the comparison stars, by summing the counts in a 5 pixel (1'.'85) 
aperture. Background is the median level in an annulus from 10- 
20 pixel (3'.'69 - 7'.'38) radius. Since VI 1 18 Ori is in a region of 
complex nebulosity, we also tried modelling the background by 
taking radial cuts through the image and extrapolating the level 
at the position of the target. This affected the magnitudes at less 
than 5% at B, and even less at longer wavelengths. 

In the near-infrared, we use 6 comparison stars in the field 
(again excluding V372 Ori). We use a 5 pixel (1'.'37) aperture. 
Background is the median level in an annulus from 10-20 pixel 
(2'.'74 - 5'.'48) radius. 

We construct a single comparison, the "superstar" from a 
weighted sum of all the comparison stars. The differential mag- 
nitude is the difference between the instrumental magnitudes of 
the target and the "superstar". 

As it happens, most of the stars in this field are known or sus- 
pected variables. Fortunately, none are as variable as VI 118 Ori. 
We assume that, by summing all the comparisons, the "super- 
star" will be less variable than any single comparison. We have 
constructed light curves using subsets of the comparison stars, 
and do not see any gross effects that can be attributed to varia- 
tions of the comparisons. Nonetheless we have not incorporated 
those stars with measured standard deviations after absolute cal- 
ibration (see below) >0. 10 mag at V into the optical "superstar". 
We caution that there is no certainty that the comparison is stable 
at the few percent level. 

A.2. Optical Absolute Calibration 

Although the SMARTS operations are not optimized for abso- 
lute photometry, single observations are made of Landolt stan- 
dard fields on each night judged photometric to establish the 
photometric zero point. These can be used to determine apparent 
magnitudes with an accuracy of better than 10%. 

On each night on which we have both observations of the 
target and of a photometric standard field, we compute the zero- 
point offset (apparent magnitude - instrumental magnitude). We 
assume the standard atmospheric extinction correction for Cerro 
Tololo. Since the standard field is observed only once per night, 
we cannot account for changes in transparency through the night. 



We do not solve for reddening terms, so the photometric solution 
consists of only the zero-point. 

We use the zero-point correction to determine the apparent 
magnitudes of all the comparison stars in the field. We could do 
this on 71 nights. In doing so, we confirm that all the comparison 
stars are variable (Table IaTI i. We co mpared these magnitudes 
to the V, Ic magnitudes published bv iHillenbrandl ([l997). The 
median offsets are 0.08 and 0.23 mag at V and Ic, respectively. 
Given that all the stars are variable members of the Orion Nebula 
population, this is acceptable agreement. 

We apply this absolute calibration on all nights (includ- 
ing non-photometric nights) to determine the apparent magni- 
tudes of the target. There are 14 nights where we overlap with 
the Villanova photometry. On these nights there are system- 
atic offsets, in the sense (SMARTS-Viflanova), of 0.16±0.10, 
0.24±0.04, and 0.34±0.04 mag in V, Rc, and Ic, respectively. 
The cause for the systematic offsets may be the lack of a color 
term in the photometric solution, exacerbated by the fact that 
most of our comparisons have V ~ Ic colors >2 mag. To allow 
comparison between the different telescope data, we have ap- 
plied the above systematic offsets to the SMARTS magnitudes 
and present the nightly average data in Table I A. 2] (available on- 
line only). 

A.3. Near-infrared Absolute Calibration 

For the JHK absolute calibration, we compare the instrumental 
magnitudes directly to the 2MASS magnitudes of the 6 com- 
parison stars. To check for target variability, we have compared 
thes e magnitudes to tho se published i n the DENIS c atalog, and 
by Ali & Depovl (1 1995b and Hillenbrand et al.l (IT998,) . Two of 
the six comparisons (2MASS J05344159-0534249 and 2MASS 
J05344219-0533036 appear variable, with variances >0.2 mag 
at K (4 observations), and 2 others (JW 94 and JW 100) appear 
variable in the optical (Table lATTT l. 

Our ability to absolutely calibrate the photometry is ulti- 
mately limited by the stability of the zero-point. We trust in 
regression to the mean, that by using at the possible compar- 
isons, on average the zero-point will be more stable that de- 
termined from any single comparison. We do not include the 
systematic uncertainty in the zero-point in our error budgets. 
Overall, the JH K flu xes compare well to those published by 
iLorenzetti et al. I (120071) . However, we have added a systematic 
correction of +0.1 mag for J and +0.2 mag for K to match 
their photometry. VI 118 Ori is intrinsically variable, as are the 
comparison stars. The only reason for applying the systematic 
offsets to our data are to simplify comparison graphically with 
the Lorenzetti et al. data. We note that the offsets do not change 
significantly the results of this paper, they merely shift the near- 
infrared data points by 0.1 to 0.2 mag, e.g., in Figures [TO] and 

m 
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Table A.l. Comparison stars for the SMARTS photometry 



Comparison B 



IV Ori 16.19±0.35 

CSV 6218 17.40±0.13 

V759 0ri 16.81±0.12 

NGC 1977 311 14.32±0.11 

VI 175 Ori 16.59±0.10 

JW94 20.82±0.35 

JW 100 20.25±0.32 

Par 1606 17.18±0.10 

CXORRS J053447.9-053544 22.00±0.99 

NSV2246 11.32±0.26 

JW33 18.75±0.32 



V 






R 






I 




5V 


51 


14.94±0. 


,26 


13 


.97±0. 


22 


13 


.19±0. 


13 






15.90±0. 


13 


14, 


.74±0. 


,13 


13 


.68±0. 


15 


0.20 


0.36 




no 


it 


^9+0 
. jzinu. 


1 n 






1 ^ 

1 D 


n if) 

U. iU 




13.24±0. 


09 


12, 


.46±0, 


,09 


11 


.73±0. 


12 






15.19±0. 


,07 


14, 


.19±0, 


,09 


13, 


.26±0. 


12 


0.08 


0.14 


19.10±0. 


17 


17, 


.37±0, 


,10 


15, 


.45±0. 


13 


0.10 


0.29 


18.48±0. 


15 


17, 


.09±0, 


,10 


15, 


.44±0. 


12 


0.03 


0.23 


15.61±0. 


,07 


14, 


.38±0, 


,08 


13, 


.22±0. 


12 


0.03 


0.12 


19.79±0. 


,32 


17, 


.76±0. 


13 


15, 


.79±0. 


12 


0.67 


0.23 


10.30±0. 


17 


9, 


.78±0, 


16 


9 


.25±0. 


13 


0.08 


0.06 


17.30±0. 


10 


16, 


.16±0, 


,09 


14, 


.84±0. 


12 


-0.03 


0.15 



Notes: 5V and 51 are the difference between these V , I determinations and those published by Hillenbrand (1997) . 
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Table A.2. SMARTS average nightly magnitudes. 



MJD 


U 


B 


V 


R 


I 


J 


H 


K 



12.74 ±0.07 
12.72 ± 0.06 



11.68 ±0.03 
11.38 ± 0.03 



18.17 ±0.05 
18.67 ±0.11 
18.41 ± 0.08 

18.32 ±0.07 
17.70 ± 0.03 



18.37 ±0.03 17.45 ±0.02 16.32 ±0.01 14.93 ± 0.01 



18.41 ± 0.03 

18.35 ± 0.03 

17.91 ± 0.02 
15.62 ±0.01 
15.68 ±0.01 
15.61 ±0.01 

15.92 ± 0.05 
15.25 ± 0.01 

14.36 ± 0.02 

14.04 ± 0.02 

14.32 ± 0.02 

14.33 ±0.01 
14.51 ±0.01 
14.44 ± 0.01 

14.66 ± 0.02 
14.31 ± 0.05 

14.73 ± 0.02 
14.24 ± 0.01 
14.21 ± 0.02 
14.31 ±0.02 
14.30 ±0.02 

14.24 ±0.02 
14.17 ±0.03 

14.09 ±0.03 
14.08 ± 0.03 

14.25 ± 0.03 
14.24 ± 0.03 

13.77 ± 0.04 

13.78 ±0.03 
13.85 ±0.11 

14.10 ±0.05 
13.90 ±0.05 
13.85 ±0.05 
13.84 ± 0.04 
13.55 ± 0.06 

14.05 ± 0.03 
13.83 ± 0.02 
13.89 ± 0.02 

13.74 ± 0.03 
13.60 ±0.02 
13.83 ±0.03 
13.95 ± 0.02 
13.95 ±0.03 
13.99 ±0.02 

13.67 ± 0.02 
13.81 ± 0.02 
13.95 ± 0.02 



17.61 ± 0.02 

17.47 ±0.02 

17.06 ±0.01 
14.70 ± 0.01 

14.78 ±0.01 
14.70 ±0.01 

14.94 ± 0.02 
14.30 ± 0.01 

13.48 ± 0.01 
13.09 ±0.01 

13.41 ±0.01 
13.40 ± 0.01 
13.58 ±0.01 
13.52 ± 0.01 
13.75 ± 0.01 

13.42 ±0.02 
13.75 ± 0.02 
13.29 ±0.01 
13.35 ± 0.02 
13.38 ± 0.02 

13.32 ± 0.02 
13.38 ±0.02 
13.47 ±0.01 
13.25 ± 0.01 
13.25 ±0.01 

13.42 ±0.01 
13.04 ± 0.02 
12.88 ± 0.03 

13.07 ±0.05 

13.33 ±0.02 
13.14 ±0.02 
13.13 ± 0.02 
13.22 ±0.01 

12.95 ± 0.03 
13.29 ±0.01 
13.06 ± 0.02 
12.97 ±0.01 
13.01 ±0.01 
12.84 ± 0.01 
12.90 ± 0.02 

13.03 ± 0.02 
13.01 ±0.03 

13.04 ± 0.02 

12.79 ±0.02 
12.90 ± 0.02 



16.32 ±0.02 

16.18 ±0.01 
15.98 ±0.01 

13.78 ±0.01 
13.87 ±0.02 

13.79 ±0.01 
13.97 ±0.02 
13.40 ± 0.01 
12.71 ±0.01 
12.35 ±0.01 

12.59 ± 0.01 

12.60 ± 0.01 
12.76 ±0.01 

12.70 ±0.01 
12.95 ±0.01 
12.60 ± 0.02 
12.94 ± 0.01 

12.51 ± 0.01 

12.59 ± 0.01 

12.60 ±0.01 
12.50 ±0.01 

12.61 ±0.01 

12.71 ±0.01 

12.52 ±0.01 
12.49 ± 0.01 
12.65 ± 0.01 
12.64 ± 0.01 
12.31 ±0.01 
12.15 ±0.01 
12.27 ±0.04 
12.54 ±0.01 
12.39 ±0.01 
12.39 ±0.01 
12.48 ± 0.01 
12.24 ± 0.02 
12.54 ± 0.01 
12.31 ±0.01 
12.24 ±0.01 
12.27 ±0.01 

12.11 ±0.01 
12.15 ±0.01 

12.26 ±0.01 
12.21 ±0.01 

12.27 ±0.01 
12.03 ±0.01 

12.12 ± 0.01 



14.91 ± 0.01 

14.86 ±0.01 
14.72 ±0.01 

12.89 ±0.01 
12.93 ±0.01 

12.87 ±0.01 
13.09 ± 0.03 
12.56 ± 0.01 

11.90 ±0.01 

11.59 ±0.01 

11.80 ±0.01 

11.81 ±0.01 
11.96 ±0.01 

11.91 ±0.01 
12.07 ±0.01 
11.90 ± 0.02 
12.12 ±0.01 

11.74 ±0.01 
11.84 ±0.01 
11.79 ±0.01 
11.69 ± 0.01 
11.81 ± 0.01 
11.95 ±0.01 

11.77 ±0.01 

11.75 ±0.01 
11.87 ±0.01 
11.89 ±0.01 

11.60 ±0.01 
11.42 ± 0.01 
11.58 ±0.03 

11.76 ±0.01 
11.69 ± 0.01 
11.67 ±0.01 
11.75 ±0.01 
11.55 ±0.01 

11.78 ±0.01 
11.60 ±0.01 
11.55 ±0.01 
11.54 ±0.01 
11.41 ± 0.01 
11.41 ± 0.01 
11.49 ±0.01 
11.45 ±0.01 
11.51 ± 0.01 
11.31 ±0.01 
11.38 ±0.01 



11.15 ± 0.04 
11.06 ±0.03 

11.32 ±0.03 



12.54 ±0.07 11.23 ±0.03 11.02 ±0.04 



12.23 ± 0.05 
12.27 ±0.05 
11.81 ±0.05 

11.71 ±0.04 
11.18 ±0.03 
11.23 ±0.03 

11.07 ±0.03 
11.41 ± 0.03 
11.23 ± 0.03 
11.39 ±0.03 
11.81 ± 0.04 
11.25 ±0.05 

11.72 ± 0.04 
11.14 ±0.03 
11.17 ±0.04 
11.21 ±0.03 
11.23 ±0.03 
11.06 ±0.03 



11.12 ±0.03 
11.17 ±0.03 
11.12 ±0.03 

10.86 ± 0.02 
10.43 ± 0.02 
10.26 ± 0.02 
10.30 ± 0.02 
10.41 ±0.02 
10.41 ±0.02 
10.58 ±0.02 
10.60 ±0.02 
10.17 ±0.03 
10.64 ± 0.02 
10.32 ± 0.02 
10.41 ±0.02 
10.39 ± 0.02 
10.34 ± 0.02 
10.38 ± 0.02 



10.56 ±0.03 
10.64 ± 0.03 
10.53 ± 0.03 

10.35 ± 0.03 
9.85 ± 0.02 
10.0 ± 0.02 

9.77 ± 0.02 
9.94 ± 0.02 
9.97 ±0.02 

10.01 ± 0.02 
9.90 ± 0.03 

9.78 ± 0.03 
10.09 ± 0.03 

9.87 ±0.02 
9.90 ± 0.03 
10.01 ±0.03 

9.92 ± 0.03 

9.93 ± 0.03 



10.92 ± 0.03 
10.69 ± 0.03 



10.21 ± 0.02 
10.04 ± 0.02 



9.68 ± 0.03 
9.58 ± 0.03 



11.18 ±0.03 
11.23 ± 0.03 



10.74 ±0.03 
10.85 ±0.03 
10.67 ±0.03 
10.60 ±0.03 
10.60 ± 0.03 
10.80 ± 0.03 



10.23 ±0.01 
10.20 ±0.01 



9.83 ±0.01 
10.04 ± 0.01 
9.98 ±0.01 
9.93 ±0.01 
9.81 ±0.01 
9.81 ±0.01 
9.87 ±0.01 



9.67 ± 0.03 
9.80 ± 0.02 



9.38 ±0.02 

9.39 ±0.02 
9.30 ±0.02 
9.43 ± 0.02 

9.29 ± 0.02 

9.30 ± 0.02 



53470 




13.00 ± 0.02 


12.19 ±0.01 


11 


46 ± 0.01 


10.90 ± 0.03 




9.20 ± 0.02 


53476 


14.17 ±0.03 


13.23 ± 0.02 


12.48 ±0.01 


11 


69 ±0.01 


10.87 ±0.03 


10.08 ±0.02 


9.50 ± 0.03 


53477 


14.22 ±0.03 


13.25 ± 0.02 


12.48 ±0.01 


11 


67 ±0.01 


10.81 ±0.03 


9.97 ±0.02 


9.51 ±0.02 


53488 


13.86 ± 0.03 


12.94 ± 0.02 


12.18 ±0.01 


11 


42 ± 0.01 


10.45 ± 0.03 


9.84 ± 0.01 


9.42 ±0.02 


53498 


13.94 ±0.14 


13.22 ± 0.04 


12.25 ±0.02 


11 


46 ± 0.01 


10.77 ±0.03 


10.25 ±0.01 


9.38 ± 0.02 


53502 




13.08 ± 0.06 












9.11 ±0.02 


53584 


14.28 ± 0.17 


13.22 ±0.01 


12.39 ±0.01 


11 


61 ±0.01 


10.82 ± 0.08 


9.98 ± 0.02 




53608 


14.46 ± 0.01 


13.49 ± 0.01 


12.67 ±0.01 


11 


86 ± 0.01 


10.93 ± 0.05 


10.13 ± 0.02 


9.57 ± 0.03 


53618 


14.43 ± 0.01 


13.44 ± 0.01 


12.62 ± 0.01 


11 


78 ± 0.01 


11.01 ± 0.06 


10.21 ± 0.02 


9.57 ±0.03 


53625 


14.89 ± 0.01 


13.94 ± 0.01 


13.09 ± 0.01 


12 


21 ±0.01 


11.44 ± 0.04 


10.44 ± 0.02 


9.83 ± 0.03 


53632 


14.50 ± 0.02 


13.51 ±0.01 


12.72 ± 0.01 


11 


86 ± 0.01 


11.05 ± 0.05 


10.20 ± 0.02 


9.51 ±0.03 


53637 


14.62 ± 0.02 


13.54 ± 0.04 










10.24 ± 0.03 


9.30 ±0.07 


53640 


14.06 ±0.01 


13.09 ±0.01 


12.31 ±0.01 


11 


52 ± 0.01 


11.01 ±0.04 


10.02 ± 0.02 


9.46 ± 0.03 


53650 


14.03 ± 0.01 


13.12 ±0.01 


12.31 ±0.01 


11 


56 ± 0.01 


10.74 ± 0.04 


10.04 ± 0.02 


9.43 ± 0.03 


53653 


14.34 ± 0.01 


13.43 ± 0.01 


12.62 ± 0.01 


11 


77 ± 0.01 


11.39 ± 0.04 


10.48 ± 0.02 


10.02 ± 0.03 



Table A.2. continued. 
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MJD 


U 


B 


V 


R 


I 


J 


H 


K 



14.15 ±0.01 

14.67 ±0.01 
14.76 ±0.01 
14.35 ± 0.01 
14.65 ± 0.01 
14.45 ± 0.01 
14.54 ± 0.01 

14.63 ±0.01 

14.96 ±0.01 

15.35 ±0.01 

15.64 ±0.01 

15.50 ±0.01 
15.43 ± 0.01 
15.69 ± 0.01 
15.54 ± 0.01 

15.36 ± 0.01 
15.08 ±0.02 

15.54 ±0.02 
15.36 ±0.02 
15.25 ±0.01 

14.97 ±0.01 

14.96 ± 0.01 
14.84 ± 0.01 
15.00 ± 0.01 

15.00 ± 0.01 

16.01 ±0.02 
16.08 ±0.02 
16.20 ±0.02 

16.16 ±0.02 
16.10 ±0.01 

15.51 ±0.01 

15.68 ±0.01 

15.17 ±0.01 

14.98 ± 0.01 
15.48 ±0.01 
15.74 ±0.01 

15.81 ±0.01 
16.00 ±0.01 
15.73 ±0.02 
15.88 ± 0.02 

15.71 ± 0.01 

15.72 ± 0.01 
15.83 ± 0.01 
15.98 ±0.01 
16.16 ±0.01 

16.23 ±0.01 

16.48 ±0.01 
16.29 ±0.01 

16.45 ± 0.02 
16.12 ± 0.02 

16.46 ± 0.02 
16.60 ± 0.01 

16.49 ± 0.01 
17.10 ± 0.02 

16.67 ±0.02 

16.82 ±0.03 

16.97 ±0.02 
16.88 ± 0.03 
16.97 ± 0.02 
17.38 ± 0.03 
17.41 ± 0.02 
17.29 ±0.05 

17.55 ± 0.03 

17.69 ± 0.15 



13.25 ±0.01 

13.73 ±0.01 
13.82 ± 0.01 
13.38 ±0.01 
13.70 ± 0.01 

13.52 ±0.01 
13.63 ± 0.01 

13.70 ± 0.01 
13.98 ±0.01 
14.51 ± 0.01 

14.74 ± 0.01 
14.63 ± 0.01 

14.53 ±0.01 

14.77 ± 0.01 

14.61 ± 0.01 
14.44 ± 0.01 

14.14 ± 0.01 

14.63 ± 0.01 
14.41 ± 0.01 
14.30 ±0.01 
14.03 ± 0.01 
14.03 ± 0.01 
13.91 ±0.01 
14.07 ± 0.01 
14.03 ± 0.01 
15.16 ± 0.02 
15.21 ± 0.02 
15.36 ± 0.02 

15.26 ± 0.01 
15.20 ± 0.01 

14.59 ±0.01 
14.73 ± 0.01 
14.28 ± 0.01 
14.07 ± 0.01 

14.64 ± 0.01 
14.89 ± 0.01 
14.87 ±0.01 
15.11 ±0.01 
14.82 ±0.01 
15.00 ±0.01 
14.81 ± 0.01 

14.85 ± 0.01 
14.95 ± 0.01 

15.15 ±0.01 

15.27 ±0.01 

15.38 ±0.01 

15.60 ±0.01 
15.47 ±0.01 
15.57 ±0.02 
15.33 ± 0.02 
15.64 ± 0.02 

15.78 ± 0.02 

15.62 ± 0.02 
16.30 ± 0.02 

15.86 ±0.02 
16.07 ±0.02 
16.13 ± 0.02 
16.06 ± 0.02 
16.09 ± 0.02 
16.53 ± 0.03 
16.55 ± 0.03 

16.68 ± 0.05 

16.69 ± 0.03 

16.71 ±0.09 



12.45 ± 0.01 
12.88 ± 0.01 
13.01 ±0.01 
12.60 ±0.01 

12.86 ±0.01 
12.69 ± 0.01 
12.81 ±0.01 

12.87 ±0.01 
13.13 ±0.01 

13.81 ±0.01 
13.73 ±0.01 
13.66 ± 0.01 
13.85 ± 0.01 

13.68 ± 0.01 
13.54 ± 0.01 
13.26 ±0.01 

13.69 ±0.01 
13.49 ±0.01 
13.42 ±0.01 

13.20 ±0.01 

13.16 ±0.01 
13.06 ± 0.01 
13.23 ± 0.01 

13.17 ±0.01 
14.26 ±0.01 
14.32 ±0.01 
14.42 ±0.01 
14.31 ±0.01 
14.26 ±0.01 
13.68 ± 0.01 
13.79 ± 0.01 
13.40 ± 0.01 

13.21 ±0.01 

13.73 ±0.01 

13.94 ±0.01 

13.95 ±0.01 
14.17 ±0.01 
13.90 ±0.01 
14.01 ±0.01 
13.84 ± 0.01 
13.95 ± 0.01 

14.00 ± 0.01 
14.19 ±0.01 
14.31 ±0.01 

14.40 ±0.01 
14.59 ±0.01 
14.49 ±0.01 

14.62 ±0.01 
14.39 ± 0.01 
14.59 ± 0.01 

14.74 ± 0.01 

14.63 ± 0.01 

15.16 ±0.01 

14.82 ±0.01 
14.99 ± 0.02 
15.04 ±0.01 
14.95 ± 0.01 

15.01 ± 0.01 

15.41 ± 0.02 

15.42 ± 0.02 

15.52 ±0.03 

15.53 ±0.02 

15.75 ± 0.05 



11.65 ±0.01 
12.02 ±0.01 
12.16 ±0.01 

11.83 ±0.01 
12.02 ± 0.01 
11.87 ±0.01 
12.00 ±0.01 

12.05 ±0.01 
12.27 ±0.01 

12.75 ± 0.01 
12.89 ± 0.01 

12.84 ±0.01 
12.79 ±0.01 

12.93 ± 0.01 

12.76 ± 0.01 

12.65 ± 0.01 

12.38 ±0.01 
12.75 ±0.01 
12.61 ± 0.01 
12.54 ± 0.01 
12.37 ±0.01 
12.35 ±0.01 

12.26 ±0.01 

12.41 ±0.01 
12.32 ± 0.01 
13.35 ±0.01 

13.42 ± 0.01 
13.51 ±0.01 

13.37 ±0.01 
13.32 ±0.01 

12.77 ±0.01 

12.85 ±0.01 

12.53 ± 0.01 

12.39 ±0.01 

12.86 ± 0.01 

13.08 ±0.01 

13.06 ±0.01 

13.24 ±0.01 
12.99 ± 0.01 
13.11 ±0.01 
12.98 ±0.01 
13.05 ±0.01 

13.09 ± 0.01 

13.25 ±0.01 

13.38 ±0.01 

13.45 ±0.01 
13.61 ±0.01 

13.54 ±0.01 

13.66 ±0.01 

13.47 ±0.01 
13.57 ±0.01 
13.75 ±0.01 

13.67 ±0.01 

14.07 ± 0.02 
13.82 ±0.01 

13.94 ± 0.02 
14.00 ± 0.02 

13.95 ±0.02 
13.98 ± 0.02 

14.27 ± 0.02 
14.27 ± 0.02 
14.35 ± 0.02 
14.38 ± 0.02 

14.48 ± 0.04 



10.93 ± 0.04 
11.23 ± 0.04 
11.73 ±0.05 
11.05 ±0.04 

11.30 ±0.04 
11.05 ±0.04 
11.21 ±0.04 
11.27 ±0.04 
11.44 ±0.05 
11.75 ± 0.07 

11.84 ±0.06 
11.80 ± 0.05 
11.95 ± 0.06 

11.84 ± 0.06 

11.72 ± 0.05 

11.48 ±0.05 

11.87 ±0.05 

11.73 ±0.05 

11.69 ±0.06 

11.50 ±0.05 

11.51 ±0.05 

11.36 ±0.05 

11.49 ± 0.05 
11.54 ±0.05 
12.15 ±0.07 
12.19 ±0.07 

11.85 ±0.08 
12.19 ±0.06 

12.03 ±0.06 
11.83 ± 0.06 

11.94 ± 0.06 

11.70 ± 0.05 
11.54 ±0.05 
11.77 ±0.06 

11.88 ±0.07 
12.09 ±0.06 

12.12 ±0.07 
11.90 ±0.06 

11.97 ±0.07 
11.79 ± 0.07 

11.98 ± 0.07 

12.04 ± 0.06 
12.09 ±0.07 
12.08 ±0.07 
12.08 ±0.07 

12.35 ±0.07 

12.37 ±0.07 

12.31 ±0.07 

12.26 ±0.07 

12.13 ±0.07 
12.08 ± 0.08 
12.42 ± 0.07 

12.27 ±0.08 

12.49 ±0.07 
12.34 ±0.08 
12.29 ±0.08 
12.51 ±0.08 
12.33 ± 0.08 
12.29 ± 0.07 

12.50 ± 0.09 
12.41 ± 0.08 

12.36 ±0.08 

12.47 ± 0.08 

12.48 ± 0.08 



10.12 ±0.01 

10.36 ±0.02 

10.86 ±0.02 
10.18 ±0.02 
10.30 ± 0.02 

10.26 ± 0.02 
10.33 ± 0.02 

10.41 ±0.02 
10.53 ± 0.02 

10.95 ±0.02 
10.91 ±0.02 
10.89 ±0.02 
10.85 ± 0.02 

10.93 ± 0.02 

10.87 ± 0.02 

10.80 ± 0.02 
10.60 ± 0.02 

10.94 ± 0.02 
10.78 ± 0.02 
10.73 ± 0.02 
10.66 ± 0.02 
10.64 ± 0.02 
10.57 ± 0.02 
10.66 ± 0.02 
10.62 ± 0.02 
11.08 ±0.03 
11.11 ± 0.03 

11.20 ± 0.03 
11.04 ± 0.03 

11.02 ± 0.02 
10.68 ±0.02 

10.81 ± 0.03 
10.70 ± 0.02 
10.64 ± 0.02 
10.91 ± 0.02 

11.03 ±0.03 

11.06 ±0.03 
11.03 ±0.03 
10.91 ±0.02 

11.00 ±0.03 

11.01 ±0.03 

10.96 ± 0.03 
10.99 ± 0.03 
11.01 ±0.03 

11.07 ±0.03 

11.08 ±0.03 
11.14 ± 0.03 
11.30 ±0.03 

11.23 ±0.03 

11.24 ±0.03 
11.14 ±0.03 

11.23 ±0.03 

11.35 ±0.03 

11.25 ±0.03 

11.37 ±0.03 

11.24 ±0.03 
11.23 ±0.03 

11.36 ±0.03 

11.21 ±0.03 
11.28 ±0.03 

11.37 ±0.03 
11.30 ± 0.03 

11.27 ±0.03 
11.32 ±0.03 

11.42 ±0.03 



9.46 ± 0.02 
9.66 ± 0.03 

10.14 ± 0.03 

9.47 ± 0.03 
9.74 ± 0.03 
9.57 ±0.03 
9.63 ± 0.03 
9.69 ± 0.03 
9.81 ± 0.03 

10.24 ± 0.03 

10.08 ± 0.04 
9.97 ± 0.03 

10.12 ±0.04 

10.08 ± 0.04 
9.95 ± 0.03 
9.83 ± 0.03 

10.22 ± 0.03 

10.06 ±0.03 

10.07 ±0.04 
9.93 ±0.03 

9.89 ± 0.03 
9.83 ± 0.03 

9.90 ± 0.03 

9.93 ± 0.03 
10.36 ±0.04 
10.34 ±0.04 
10.10 ±0.04 
10.22 ±0.04 

10.20 ± 0.04 
9.79 ± 0.03 

9.94 ± 0.04 
9.81 ± 0.03 
9.89 ± 0.03 

10.18 ±0.04 
10.07 ±0.05 
10.17 ±0.04 

10.21 ±0.04 
10.20 ± 0.04 
10.33 ± 0.04 

10.30 ± 0.04 

10.31 ± 0.04 
10.24 ± 0.04 
10.36 ± 0.04 
10.41 ±0.04 
10.44 ± 0.04 

10.64 ± 0.04 
10.63 ± 0.04 

10.58 ± 0.04 
10.46 ± 0.05 
10.57 ± 0.05 

10.65 ± 0.05 
10.78 ± 0.05 

10.66 ± 0.05 

10.80 ± 0.05 

10.83 ±0.05 
10.69 ±0.05 
10.69 ±0.05 

10.59 ± 0.05 
10.68 ± 0.05 

10.81 ± 0.05 
10.78 ± 0.05 
10.62 ± 0.06 
10.77 ± 0.05 
10.85 ± 0.05 



Table A.2. continued. 
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MJD 


U 


B 


V 


R 


I 


J 


H 


K 



17.80 ±0.05 
18.50 ± 0.29 
18.29 ±0.29 

18.56 ±0.04 
16.50 ±0.01 
16.85 ±0.01 

16.75 ±0.01 

16.79 ±0.01 

16.57 ±0.01 
16.66 ± 0.01 
16.54 ± 0.01 
16.78 ±0.01 

16.80 ±0.01 
16.94 ± 0.01 
16.94 ± 0.02 
17.08 ± 0.02 

16.76 ± 0.01 
16.80 ±0.01 
17.71 ± 0.04 
17.82 ± 0.03 
18.42 ± 0.04 
lcS.35 ± o.o:i 



16.98 ±0.07 
17.23 ±0.13 

17.85 ± 0.28 

17.58 ± 0.03 
15.61 ±0.01 

15.99 ± 0.01 

15.91 ±0.01 

15.92 ±0.01 
15.71 ±0.01 
15.79 ±0.01 
15.67 ±0.01 
15.92 ± 0.01 
16.12 ±0.01 
16.11 ±0.01 
16.10 ±0.01 
16.19 ±0.02 
15.91 ±0.01 
15.98 ±0.01 

16.86 ±0.03 
16.95 ±0.02 
17.55 ± 0.03 
17.48 ± 0.02 



15.76 ± 0.04 
16.33 ± 0.08 
16.22 ± 0.09 

16.26 ± 0.03 
14.67 ± 0.01 
14.93 ±0.01 
14.84 ± 0.01 
14.86 ± 0.01 
14.69 ± 0.01 
14.84 ± 0.01 
14.71 ±0.01 

14.91 ±0.01 

15.11 ±0.01 
15.16 ±0.01 

15.12 ±0.01 
15.18 ± 0.02 

14.92 ± 0.01 
14.97 ±0.01 
15.71 ±0.03 
15.75 ± 0.02 
16.18 ± 0.03 
10.18 ± 0.02 



14.53 ±0.04 
14.79 ±0.05 

14.84 ±0.05 

14.89 ± 0.08 
14.81 ± 0.02 
13.74 ± 0.01 
13.95 ± 0.01 

13.85 ±0.01 

13.90 ±0.01 
13.74 ±0.01 

13.86 ±0.01 
13.76 ±0.01 
13.95 ± 0.01 
14.08 ± 0.01 
14.15 ± 0.01 
14.08 ± 0.01 
14.15 ±0.01 
13.90 ±0.01 
13.98 ±0.01 
14.52 ±0.01 
14.55 ±0.01 
14.81 ± 0.02 
14.79 ± 0.01 



12.42 ± 0.13 
12.60 ± 0.10 
12.70 ± 0.08 
12.73 ±0.18 
12.78 ±0.10 
12.42 ± 0.08 

12.48 ± 0.09 
12.45 ± 0.08 
12.51 ± 0.09 

12.35 ±0.07 

12.36 ±0.07 
12.33 ±0.07 
12.31 ±0.08 
11.75 ±0.06 
12.47 ± 0.07 
12.45 ±0.07 
12.51 ± 0.09 
12.38 ±0.07 

12.49 ± 0.08 
12.66 ± 0.09 
12.62 ± 0.08 
12.69 ± 0.09 
12.05 ± 0.08 



11.41 ±0.04 
11.60 ±0.04 
11.60 ±0.03 
11.83 ±0.12 
11.60 ± 0.04 

11.41 ± 0.03 
11.45 ± 0.03 

11.40 ±0.03 

11.42 ±0.03 
11.33 ±0.02 
11.31 ±0.03 
11.30 ±0.02 

11.38 ± 0.03 
11.33 ±0.03 

11.43 ± 0.02 
11.37 ±0.03 
11.43 ±0.03 

11.39 ±0.03 

11.41 ±0.03 

11.51 ±0.04 

11.52 ±0.03 
11.62 ± 0.03 
11.00 ± 0.03 



10.52 ± 0.07 
11.06 ±0.06 

10.82 ± 0.03 
10.91 ±0.17 
11.16 ±0.06 
10.96 ± 0.05 
10.98 ± 0.06 
10.93 ± 0.05 
10.96 ±0.05 

10.77 ±0.04 
10.81 ± 0.04 

10.78 ± 0.04 
10.81 ±0.05 
10.78 ± 0.04 

10.83 ± 0.04 
10.87 ± 0.04 
10.80 ± 0.05 
10.77 ± 0.04 
10.80 ± 0.04 
11.08 ±0.05 
11.06 ±0.05 
11.21 ± 0.05 
11.19 ± 0.05 
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Table A.3. Villanova magnitudes. 



MJD V R I 



53398 


13. 


,66 


± 


0, 


.05 


12, 


.81 


± 


0. 


,06 


12. 


,00 


± 


0. 


,02 


53399 


13, 


,55 


± 


0, 


,02 


12, 


.71 


± 


0, 


,03 


11, 


,91 


± 


0, 


,01 


53401 


13, 


,62 


± 


0, 


.07 


12, 


.75 


± 


0, 


,05 


11, 


,95 


± 


0, 


,04 


53402 


13, 


,70 


± 


0, 


.06 


12, 


.81 


± 


0, 


,05 


12, 


,02 


± 


0, 


,01 


53403 


13, 


,81 


± 


0, 


.06 


12, 


.92 


± 


0, 


,05 


12, 


,11 


± 


0, 


,02 


53406 


13, 


,64 


± 


0, 


.03 


12, 


.78 


± 


0, 


,05 


11, 


,97 


± 


0, 


,02 


53418 


13, 


,58 


± 


0, 


.33 


12, 


.49 


± 


0, 


.09 


11, 


.71 


± 


0, 


.04 


53428 


13, 


,09 


± 


0, 


.07 


12, 


.28 


± 


0, 


.04 


11, 


.53 


± 


0, 


.03 


53435 


12, 


,77 


± 


0, 


.02 


11, 


.95 


± 


0, 


.01 


11, 


.18 


± 


0, 


.05 


53436 


12, 


,57 


± 


0, 


.02 


11, 


.78 


± 


0, 


.02 


11, 


.05 


± 


0, 


.01 


53444 


13, 


,23 


± 


0, 


.03 


12, 


.37 


± 


0, 


,03 


11, 


,57 


± 


0, 


,03 


53445 


13, 


,21 


± 


0, 


.03 


12, 


.36 


± 


0, 


,03 


11, 


,56 


± 


0, 


,03 


53621 


13, 


,58 


± 


0, 


.04 


12, 


.68 


± 


0, 


,03 


11, 


,86 


± 


0, 


,02 


53625 


13, 


,99 


± 


0, 


.03 


13, 


.07 


± 


0, 


,05 


12, 


,20 


± 


0, 


,02 


53634 


13, 


,52 


± 


0, 


.06 


12, 


.64 


± 


0, 


,03 


11, 


,80 


± 


0, 


,04 


53641 


13, 


,10 


± 


0, 


.04 


12, 


.26 


± 


0, 


.02 


11, 


.50 


± 


0, 


.02 


53659 


13, 


,49 


± 


0, 


.04 


12, 


.62 


± 


0, 


.03 


11, 


.79 


± 


0, 


.03 


53674 


13, 


,44 


± 


0, 


.02 


12, 


.60 


± 


0, 


.02 


11, 


.82 


± 


0, 


.02 


53707 


14, 


,09 


± 


0, 


.03 


13, 


.17 


± 


0, 


.04 


12, 


.34 


± 


0, 


.03 


53711 


14, 


,82 


± 


0, 


.09 


13, 


.89 


± 


0, 


,08 


12, 


,99 


± 


0, 


,03 


53712 


14, 


,74 


± 


0, 


.13 


13, 


.81 


± 


0, 


,06 


12, 


,89 


± 


0, 


,04 


53718 


14, 


,45 


± 


0, 


.12 


13, 


.57 


± 


0, 


,11 


12, 


,69 


± 


0, 


,05 


53748 


15, 


,12 


± 


0, 


.22 


14, 


.20 


± 


0, 


,16 


13, 


,28 


± 


0, 


,11 


53751 


15, 


,12 


± 


0, 


.34 


14, 


.36 


± 


0, 


,19 


13, 


,41 


± 


0, 


,09 


53754 


14, 


,49 


± 


0, 


.16 


13, 


.55 


± 


0, 


.04 


12, 


.68 


± 


0, 


.04 


53755 


14, 


,61 


± 


0, 


.04 


13, 


.69 


± 


0, 


.06 


12, 


.78 


± 


0, 


.03 


53762 


14, 


,16 


± 


0, 


.08 


13, 


.30 


± 


0, 


.08 


12, 


.44 


± 


0, 


.06 


53773 


14, 


,99 


± 


0, 


.08 


14, 


.11 


± 


0, 


.07 


13, 


.19 


± 


0, 


.05 


53774 


14, 


,78 


± 


0, 


.11 


13, 


.92 


± 


0, 


,07 


13, 


,00 


± 


0, 


,02 


53785 


14, 


,83 


± 


0, 


.12 


13, 


.94 


± 


0, 


,08 


13, 


,06 


± 


0, 


,06 


53790 


15, 


,03 


± 


0, 


.07 


14, 


.12 


± 


0, 


,06 


13, 


,18 


± 


0, 


,03 


53793 


15, 


,39 


± 


0, 


.14 


14, 


.48 


± 


0, 


,10 


13, 


,52 


± 


0, 


,03 


53795 


14, 


,87 


± 


0, 


.08 


14, 


.09 


± 


0, 


,05 


13, 


,19 


± 


0, 


,03 


53800 


15, 


,23 


± 


0, 


.10 


14, 


.37 


± 


0, 


.08 


13, 


.42 


± 


0, 


.05 


53802 


15, 


,16 


± 


0, 


.08 


14, 


.24 


± 


0, 


.05 


13, 


.36 


± 


0, 


.03 
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Appendix B: Details on thie Spitzer data 
reduction 

B. 1. IRAC 

We used the MOsaicker and Point source Extractor (MOPEX) 
software release of June 2007. We started from the Basic 
Calibrated Data (BCD) individual images to produce mosaics 
with a (native) pixel size of 1'.'22. For the data taken in sub-array 
mode (PIDs 3716 and 41019), we collapsed the 3-dimensional 
image BCDs into 2-dimensional BCDs by using, for each pixel, 
the median of the 64 planes. For the uncertainty BCD files, we 
used the standard deviation of the 64 planes of the BCD file in- 
stead of using the input uncertainties. This method allows the 
removal of most particle hits and the ingestion of the collapsed 
BCDs in MOPEX. We then used aperture photometry, centered 
on VI 118 Ori, using an extraction radius of 10 pixels (12'.'2), 
and a concentric annulus of radii 12 and 15 pixels (we used radii 
of 10 and 12.5 pixels for the IRAC4 band at 8.0 ^m due to the 
strength of the background intensity at this wavelength for radii 
larger than 12.5 pixels. Indeed, the background is dominated by 
emission due to the nearby Herbig Ae star V372 Ori; this effect 
is less prominent in the other IRAC bands, but the presence of 
the Herbig star limited us to outer radii of 15 pixels). No aperture 
correction is needed for an extraction circle radius of 10 pixels 
(according to the IRAC data handbook v3.0). 

B.2. MIPS 

We used the BCD files where V 1 1 1 8 Ori was on the detector 
as input files for the MOPEX pipeline. A native pixel size of 
2'.'5 was used to create the mosaic, and we used an extraction 
radius of 13" centered on VI 1 18 Ori to derive the aperture pho- 
tometry. We used an annulus of radii 15" and 20" for the back- 
ground, and finally used an aperture correction factor of 1.17 
(from Table 3.12 of the MIPS data handbook, v. 3.3.1). We also 
investigated the background with a different method, i.e., to cal- 
culate the background in a nearby region using a circle with ra- 
dius 13". Indeed, the MIPS24 background near VI 118 Ori is 
highly inhomogeneous. Notice that, while the uncertainty in the 
flux density is of the order 0.5 mJy, we estimate the true uncer- 
tainty, mainly due to the inhomogeneity of the background and 
the difficulty to find a representative background region, to be 
closer to 15-20 mJy. 

B.3. IRS 

We extracted the IRS spectra using the Spitzer IRS Custom 
Extraction (SPICE) v2.1.2 software and the post-BCD co-added 
2-D images (pipeline version S 15.3.0 for SL/SH and S 17.2.0 for 
LH/LL for the outburst data; S18.5.0 for all post-outburst spec- 
tra). We investigated in detail the methodology to subtract the 
background in the IRS spectra. Indeed, for the low-resolution 
modules, we have investigated different techniques to estimate 
the background contamination. First, we used the standard tech- 
nique that takes the image from nod A, subtracts the nod B im- 
age, and then the net ( A-B) spectrum of V 1 1 1 8 Ori is extracted 
with SPICE. The same is then done for the net (B-A) spectrum, 
and both net spectra are averaged. Unfortunately, this technique 
is unsatisfactory due to the inhomogeneity of the background 
emission along the cross-dispersion axis of the slit. A second 
technique consisted in using the images taken for the SL:^ nods 

' For clarity, we explain the second technique for SL2 only in the 
main text. For the SLl data, simply permute SL2 and SLl in the text. 



and use the source-free region in the SLl image for the back- 
ground of the SLl nods. This method removes any contamina- 
tion by VI 118 Ori in the SLl image, but the downside is that 
the region of the sky covered by the SLl slit (during the SL2 
observations) are away from VI 118 Ori. Thus, the measured 
background may not represent the background near VI 1 18 Ori. 
A third technique consisted in using a background region near 
VI 1 18 Ori obtained during the observed nod. This technique al- 
lows to get a better estimate for the background near VI 1 18 Ori, 
but has the disadvantage that VI 118 Ori may contaminate the 
background, especially at longer wavelengths, since the stan- 
dard extraction width increases with increasing wavelength, to 
take into account the increasing size of the point spread function 
of a point source. For the February 2005 observation, technique 
#2 proved better for SLl, while technique #3 was better for SL2. 
We used the PAH emission at 6.2 and 11.3 /im (coming mostly 
from the diffuse background emission) to check that the back- 
ground emission was adequately removed. For the March 2005 
observation, we used technique #3 for both SL slits, as technique 
#2 gave an excess in the PAH at 11.3 /im in the SLl slit. For the 
post-outburst observations, we used technique #3 for both SL 
slits and also for the LL slits (in this case, the background sub- 
traction is more difficult for A > 30 mic). 

For the high-resolution modules, we used the standard ap- 
proach to average the spectra from both nods using the full slit 
aperture. For the post-outburst spectra, we used the accompa- 
nying background spectra to obtain the background-subtracted 
spectra. This procedure worked nicely for the SH module (i.e., 
up to 19.35 iim), as demonstrated by the good subtraction of 
the PAH emission at 11.3 ^m. Notice that the H2 - 5(1) 
A17.0348 fim line was apparently completely subtracted, while 
the H2 - 5(2) A12.2786 /im line is faintly detected. We 
also detect [Nell] A12.8136 fim CP1/2 - ^^3/2) an [S III] 
A18.7130 iim (^P2 - "^Pi) in the background-subtracted spec- 
trum (we will come back to this issue below). Furthermore, the 
silicate feature and continuum flux up to 14 /im are consistent 
with the SL flux, suggesting that the background subtraction 
was not too far off. For the LH module, we subtracted only 95% 
of the background flux for the post-outburst data. Indeed, at the 
wavelengths covered by the LH, the background dominates the 
emission and is highly position-dependent. Although our back- 
ground pointing was close to VI 118 Ori, there is considerable 
inhomogeneity in the diffuse emission of the Orion nebula, as 
observed in the 8.0 and 24 /tm images. Figure IB72] shows the 
MIPS 24 /tm images together with the "field-of-views" of the 
IRS slits during the first outburst observation and in late 2008. 
We used the [Si ll] A34.815 /im CP3/2 - ^^1/2) as the proxy for 
background subtraction for the LH module, and noticed that a 
scaling factor of 0.95 was better than 1.0. We detect, however, 
excess flux in the [S III] A33.481 /im (^Pi - ^Pg) line and faint 
excess in the H2 — 5(0) A28.2188 fim line (in contrast to 
the S(l) line in the SH spectrum). The SH and LH spectra rel- 
atively well at 19.4 /im, suggesting again that the background 
subtraction was adequate. 

We provide, however, a word of caution about the accuracy 
of the continuum flux level (especially above 30 /im) and in the 
reality of the detected flux excess in some lines, especially in the 
LH wavelength region. Indeed, we have scaled the background 
LH spectrum using the [Si ll] line. However, we remind that 
VI 118 Ori is located at the center of the Orion nebula, where 
strong ionization and excitation of the interstellar medium takes 
place. Strong variations of the diffuse background emission take 
place, and they may differ for the continuum emission and for the 
line emission (mostly [Si ll] , [S III], [Ne ll], and H2). Therefore, 
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using a specific emission line to check that the background sub- 
traction is adequate may result in an incorrect subtraction of the 
continuum emission of the nebula and even in the subtraction 
of another emission line! Nevertheless, the similarities of the 
continuum shape in the low-resolution and high-resolution mod- 
ules (at least below 35 fim) for the post-outburst data, the pres- 
ence of line excess in [S III] in the LL spectrum and in both SH 
and LH spectra are suggestions that the background subtraction 
was overall accurate. The presence of [Ne ll] in the SH module 
but not in the SL module demands, however, deeper analysis to 
confirm the detection of the line. We degraded the SH and LH 
spectra down to SL and LL resolution (Fig. lB.ll available online 
only). The continuum shapes are well matched. But the degraded 
SH data covering the [Ne ll] line shows a strong line which is not 
compatible with the SL data at this wavelength. This strongly 
suggests that the SH background [Ne ll] line flux was too faint 
compared to the line flux near VI 1 18 Ori and that the measured 
SH line flux is of background origin (from the Orion nebula or 
from diffuse emission close to VI 1 18 Ori). In the case of H2 at 
12.28 ^m, the degraded SH spectrum is consistent with the SL 
spectrum: the contrast between the line and the continuum is too 
faint to confirm the presence of this line in the SL spectrum. The 
same comment applies to the other molecular hydrogen line at 
28.2 /im. Thus, we cannot confirm that this molecule is detected 
in VI 118 Ori's spectrum. On the other hand the [S III] lines in 
the degraded SH and LH spectra have similar peak flux densi- 
ties as in the LL spectrum, indicating that these lines originate 
from the immediate vicinity of V1118 Ori. In summary, we be- 
lieve in the detection of the [S III] lines, while we have doubts 
for [Ne II] and the H2 lines, but we provide nevertheless the line 
fluxes from the high-resolution module spectra in Table[3] In any 
case, in view of Spitzer spatial resolution, higher spatial resolu- 
tion observations would be required to confirm the origin of the 
excess line emission. 

We have further investigated whether the post-outburst back- 
ground observations could be used for the outburst SH and LH 
data. In principle, the zodiacal light dominates the continuum 
background emission below about 40 /im, while the interstel- 
lar medium emission dominates above. The zodi contribution 
varies as a function of time. We have used the Spitzer Planning 
Observations Tool (SPOT) to determine the contribution at the 
three different epochs of the IRS observations and found that 
they were of similar level (about 20-25 MJy sr^^ from 15 to 
24 /xm). For the SH spectrum, we preferred to use the post- 
outburst spectrum without scaling factor, since the PAH feature 
and the H2 line S(2) were well canceled (the S(l) line is slightly 
oversubtracted). For the LH spectrum, we used a scaling factor 
equivalent to 0.85, i.e., about the value of the ratio of the zodi- 
acal light contribution at 24 at in February-March 2005 and 
November 2008. This ratio also cancels out relatively well the 
[Si II] line while the H2 S(0) line still remains detected. With 
this procedure, the resulting spectrum shows an increase of the 
flux with longer wavelengths, perhaps coming from an envelope, 
but this increase is not detected in the post-outburst data, cast- 
ing some doubt on the accuracy of the continuum shape in the 
SH and LH outburst spectra above 14 /tm. Since we have no 
outburst LL data to confirm this shape, we prefer to err on the 
safe side and claim that the high-resolution continuum shape for 
A > 14 /im during the outburst is unreliable. The strengths of 
the emission lines may also be affected by improper background 
reduction; indeed the flux levels are generally larger than in the 
post-outburst spectrum (except for [Nell], but in this case, the 
line is probably not originating in the direct vicinity of VI 118 



Ori). Line fluxes from the high-resolution module data in out- 
burst are also given in Table |3] 
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Fig.B.l. Spitzer IRS high-resolution post-outburst spectra (thick 
curve) degraded to the resolution of the SL and LL modules. The 
latter data are shown as well (thin curve). Notice the good agree- 
ment in the continuum and for most lines, except for [Ne ll]. 
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Fig.B.2. Spitzer MIPS 24 image taken in March 2004 with 
overlays of the IRS modules during the observations on 2005 
February 18 (top) and on 2008 November 14 after the outburst 
(bottom). Note that we only show the first position of the two 
nods, for clarity. Both nods for the background SH/LH obser- 
vations are shown in the bottom figure. A 1' scale is shown on 
both North-orientated images. The image scale is linear from 
to 250 MJy sr-i. 



